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Abstract 


The  structure  of  the  near-bottom  velocity  and  density  fields  was 
observed  for  3  months  with  a  fixed  velocity/temperature  measuring  array 
jn  the  Hatteras  Abyssal  Plain.  The  velocity  signal  was  examined  for 
structure  coherent  with  the  observed  mixed  layers.  Velocity  fluctuations 
above  1  cph  increase  in  energy  near  the  bottom,  expecially  within  bottom 
mixed  layers.  The  frequency  and  energy  of  these  fluctuations  is  con¬ 
sistent  with  the  expected  properties  of  boundary  layer  turbulence.  The 
turbulence,  as  measured  by  these  velocity  fluctuations,  is  modulated 
on  tidal  and  inertial  time  scales  and  extends  intermittently  throughout 
the  bottom  mixed  layer.  The  clpckwise  near  inertial  velocity  fluctua¬ 
tions,  presumably  due  to  internal  waves,  also  show  structure  coherent 
w.th  the  observed  mixed  layers.  Their  energy  decreases  near  the  top  of 
the  mixed  layer  with  little  phase  change.  Within  the  mixed  layer  the 
phase  begins  to  lead,  with  the  phase  increasing  downward.  The  near 
inertial  anticlockwise  velocity  fluctuations  show  far  less  coherence  with 
th<’  mixed  layer  structure;  higher  frequency  internal  wave  band  velocity 
fluctuations  show  very  little  coherence  with  the  mixed  layer  structure. 
The  characteristic  boundary  layer  velocity  signal  in  both  the  high 
frequency  and  near  inertial  bands  commonly  extends  throughout  the  mixed 
layer,  often  significantly  above  the  estimated  turbulent  Ekman  layer 
height.  These  observations  are  inconsistent  with  a  steady  turbulent 
Ekman  layer  model  of  the  boundary  layer. 
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The  interaction  of  the  internal  wave  field  and  the  benthic  boun¬ 
dary  layer  above  a  flat  bottom  is  investigated  using  a  diagnostic  model. 
Assuming  the  internal  waves  to  be  large  compared  with  the  boundary 
layer,  the  turbulent  stresses  acting  on  the  waves  and  the  work  done  by 
the  stresses  on  the  waves  can  be  computed.  The  data  indicates  obser¬ 
vable  energy  transfer  only  for  near  inertial  waves.  The  boundary  layer 
is  estimated  to  absorb  .004  -  .025  ergs/cm2 /s  from  the  near  inertial 
internal  wave  field.  This  is  far  less  than  estimated  by  Leaman  (1976) 
and  suggests  that  the  benthic  boundary  layer  on  a  flat  bottom  plays  a 
minor  role  in  dissipating  internal  wave  energy.  This  is  also  much  less 
than  the  total  energy  dissipation  in  the  boundary  layer,  suggesting  that 
the  boundary  layer  is  primarily  driven  by  low  frequency  motions.  A 
simple  slab  model  with  a  linearized  quadratic  drag  law  qualitatively 
explains  the  observed  boundary  layer  near  inertial  clockwise  velocity 
structure  and  energy  flux.  Thesis  Advisor:  Laurence  Armi. 
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Chapter  1.  Introduction 

Observations  close  to  the  ocean  bottom  show  a  distinct  region, 
the  benthic  boundary  layer,  in  which  bottom  generated  shear  turbulence 
appears  to  be  dynamically  significant.  Bowden  (1978)  reviews  previous 
work  on  this  subject.  Oceanographically,  the  benthic  boundary  layer 
may  be  significant  in  damping  oceanic  motions,  and  mixing  the  oceanic 
interior.  Fluid  dynamically,  it  is  a  useful  example  of  a  turbulent 
boundary  layer  influenced  by  rotation  and  stratification,  but  without 
significant  heating  or  cooling.  Over  large  regions  of  the  ocean  a 
bot-om  mixed  layer,  vertically  uniform  in  potential  temperature,  salinity 
and  light  scattering  (Armi  and  Millard,  1976;  Caldwell,  1976;  Weatherly 
and  Niiler,  1974)  is  present.  This  layer  is  believed  to  be  mixed  by  the 
boundary  layer  turbulence.  In  a  recent  study  Armi  and  D’Asaro  (1980) 
describe  the  three  dimensional  density  structure  of  the  benthic  boundary 
layer  over  the  Hatteras  Abyssal  Plain  in  the  western  North  Atlantic. 

The  bottom  mixed  layer  is  found  to  extend  5-60  m  above  the  abyssal  plain 
with  typical  heights  of  20-30  m.  The  layer  exhibits  considerable 
variability  on  time  scales  of  days  and  longer,  and  horizontal  space 
scales  of  tens  of  kilometers  and  larger.  Multiple  mixed  layers  are 
common.  These  are  presumably  caused  by  the  intrusion  of  the  bottom 
mixed  layer  into  the  interior.  Benthic  fronts  are  also  observed.  The 
observed  density  variations  at  a  fixed  point  appear  to  be  primarily  due 
to  advection  of  horizontal  density  structure  by  the  velocity  field. 

This  thesis  will  analyze  the  velocity  data  taken  during  the  same 
experiment.  Studies  of  the  benthic  boundary  layer  have  often  concen- 
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trated  on  the  structure  of  the  mean  flow  and  turbulence  in  the  bottom 
part  r-f  Lhe  boundary  layer.  In  this  region  the  benthic  boundary  -.yer 
appears  similar  to  a  two  dimensional  turbulent  boundary  layer  and  the 
velocity  profile  can  be  used  to  estimate  the  bottom  stress.  (Wimbush 
and  Munk,  1970;  Weatherly,  1972;  Caldwell  and  Chriss,  1079).  The 
measurements  described  here  are  not  sufficiently  close  to  the  bottom  to 
duplicate  these  calculations.  They  instead  span  the  entire  mixed  layer, 
extending  into  the  interior,  so  that  the  velocity  structures  associated 
with  the  mixed  layer  can  be  observed. 

The  turbulent  Ekman  layer  provides  a  useful  starting  point  for  the 
analysis  of  the  mixed  layer  velocity  structure.  The  unstratified  tur¬ 
bulent  bounuary  layer  driven  by  a  steady  geostrophic  flow  in  a  rotating 
coordinate  system  will  here  be  called  a  turbulent  Ekman  layer.  Its 
structure  has  been  investigated  theoretically  using  asymptotic  matching 
(Csanady,  1967;  Blakadar  and  Tennekes,  1968),  laboratory  simulations 
(Caldwell  et  ad.,  1972;  Howroyd  and  Slawson,  1975),  and  the  analysis  of 
atmospheric  and  oceanic  boundary  layers  under  neutral  conditions 
(r.cttau,  1950;  McPhee  and  Smith,  1976).  One  result  of  these  studies 
is  that  the  turbulent  Ekman  layer  extends  only  to  a  height  of  roughly 
iS  v  .4  u*/f  where  pu*2  is  the  wall  stress.  The  turbulent  stresses 
and  velocity  fluctuations  diminish  rapidly  above  this  height  (Caldwell 
et  al. ,  1972;  McPhee  and  Smith,  1976).  The  benthic  boundary  layer  is 
certainly  more  complex  than  a  turbulent  Ekman  layer.  One  of  the  main 
aim-,  of  this  thesis  is  determining  the  relevance  of  the  turbulent  Ekman 
layer  and  related  models  to  the  benthic  boundary  layer. 
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This  thesis  divides  into  two  sections  by  content.  Chapters  3-6 


describe  the  observed  mixed  layer  velocity  structures.  High  frequency 
(chapter  4)  and  near  inertial  frequency  (chapter  5)  motions,  which  show 
the  strongest  structure  are  emphasized.  These  results  are  discussed  in 
chapter  6.  Chapters  7-10  analyze  the  near  inertial  motions  using  a 
diagnostic  model  of  the  internal  wave/boundary  layer  interaction  which 
allows  the  energy  flux  between  the  waves  and  the  boundary  layer  to  be 
computed.  This  interaction  is  then  modeled  using  a  simple  slab  model 
(chapter  9)  and  the  results  discussed  (chapter  10) . 
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Chapter  2.  The  experiment  and  the  instruments 

A  3-month  time  series,  May  18,  1977  to  August  18,  1977,  of 
velocity  and  temperature  was  measured  by  a  bottom  mooring  (Figure  1)  de¬ 
ployed  on  the  Hatteras  Abyssal  Plain  near  28°N,  70°  30'W  by  the  Moored 
Array  Project  of  the  Woods  Hole  Oceanographic  Institution  (W.H.O.I.). 

The  Hatteras  Abyssal  Plain  is  extremely  flat  with  a  slope  of  only 
2  *  10  '  over  hundreds  of  kilometers  (Bush,  1976) .  The  mooring 
(W.H.O.I.  621)  contained  vector-averaging  current  meters  (VACM)  spaced 
so  as  to  span  the  bottom  mixed  layer.  These  sampled  average  velocity 
and  temperature  every  7  1/2  minutes.  A  detailed  description  of  the  ex¬ 
periment  and  the  instrument  calibration  can  be  found  in  (Spencer, 

D'Asaro  and  Armi,  198l) • 

The  VACM  instruments  used  in  this  experiment  employ  a  Savonius 
rotor  and  a  vane  to  measure  horizontal  current  (McCullough,  1975) . 

Bryden  (1976)  estimates  an  rms  VACM  direction  bias  of  2°.  The  measured 
currents,  however,  show  consistent  directional  differences  larger  than 
this.  A  more  detailed  analysis  (Spencer,  D'Asaro  and  Armi,  1981  )  reveals 
that  instruments  1  and  4  are  consistently  7°  clockwise  of  instruments  2 
and  3  from  early  June  up  to  July  19.  During  this  period  the  mooring 
rotation,  as  measured  by  the  VACM  compasses  is  minimal.  On  July  19  and 
20  the  mooring  orientation  twice  changes  abruptly.  Instrument  2,  for 
example,  undergoes  a  net  rotation  of  approximately  180°  in  less  than 
the  sampling  interval.  The  directional  differences  disappear  after 
July  20.  This  strongly  suggests  that  the  shears  resulting  from  these 
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Figure  1 


Mooring  diagram.  All  vector-averaging  current  meters  (VACM)  measure 
vector  mean  current  and  temperature. 


directional  differences  should  be  interpreted  as  instrumental  errors  of 
unknown  origin.  The  VACMs  .n  this  experiment  are  therefore  considered 
to  have  possible  direction  biases  of  up  to  7°  per  pair,  with  this  bias 
likely  to  be  variable  in  time.  These  biases  dominate  the  low  frequency 
shears  observed  ir.  this  experiment,  and  make  their  scientific  inter¬ 
pretation  difficult. 

The  measured  VACM  temperature  variations  have  a  precision  of  better 
than  a  millidegree  (±  10  5°C).  The  absolute  temperature  is,  however, 
less  accurately  measured  and  has  been  calibrated  in  situ  as  described  in 
Armi  and  D'Asaro  (1980)  and  Spencer,  D'Asaro  and  Armi  (1981).  The  cor¬ 
relation  between  salinity  and  potential  temperature  in  this  region,  as 
determined  by  both  historical  measurements  (MODE  Group,  1978)  and  ex¬ 
tensive  CTD  measurements  during  this  experiment,  is  exact  to  within  the 
salinity  measurement  error.  Thus  temperature  alone  gives  a  measurement 
of  density  at  least  as  accurate  as  temperature  and  salinity  together. 
Potential  temperature  differences  will  thus  be  interpreted  as  potential 
density  differences  throughout  this  paper  using  the  formula  gAp/p  =  .  1A0 
(Armi  and  D'Asaro,  1980). 
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Chapter  3. 


Spectral  analysis 


Figure  2  shows  the  horizontal  kinetic  energy  spectra  from  all  7 
instruments.  These  spectra  are  similar  to  those  from  other  measurements 
in  this  area  showing  peaks  at  the  inertial  and  tidal  frequencies  with 
a  red  spectrum  at  higher  frequencies.  The  normalized  energy  at  .2  cph, 

3  cm*’/s2/eph2 ,  is  within  the  range  found  by  Wunsch  (1976)  for  western 
North  Atlantic  measurements.  The  spectral  slope  at  .2  cph,  approxi¬ 
mately  -3,  is  steeper  than  the  typical  slope  of  -1.7  found  by  Wunsch 
(1976) .  Figure  3  shows  the  horizontal  kinetic  energy  spectrum  from  the 
55  m  instrument,  here  divided  into  clockwise  and  anticlockwise  com¬ 
ponents  (Caiman,  1978) .  The  clockwise  energy  is  larger  than  the  anti¬ 
clockwise  energy  between  che  inertial  and  Vaisala  frequency,  as  is 
typical  of  previous  measurements  (Fu,  1980) . 

Several  features  of  the  mean  velocity  spectra  in  figure  2  show 
distinct  changes  across  the  array.  The  near  bottom  instruments  measure 
significantly  more  horizontal  kinetic  energy  at  frequencies  above  1  cph. 
In  chapter  4  this  will  be  attributed  to  bottom  generated  turbulent 
velocity  f luctuat ions .  The  near  bottom  instruments  also  measure  sig¬ 
nificantly  less  horizontal  kinetic  energy  near  the  inertial  frequency. 

In  chapter  6  this  is  attributed  to  turbulent  stresses  within  the  mixed 
layer.  In  both  cases  the  energy  level  varies  monotonically  with  distance 
from  the  bottom.  Except  at  these  frequencies  the  mean  spectra  from 
various  instruments  do  not  significantly  differ. 
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FREQUENCY  ,  CYCLES /HRS. 


Figure  3 

Clockwise  (solid)  and  anticlockwise  (dotted)  kinetic  energy  spectra  for 
55  m  instrument.  Each  spectrum  is  the  average  of  13  Hanned  overlapped 
pieces  with  higher  frequencies  band  averaged.  The  95%  confidence  limits 
and  the  inertial  (f),  tidal,  and  approximate  Vaisala  (N)  fre¬ 

quencies-  are  shown. 
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Figure  4  shows  the  east/east  and  north/north  cross  spectra  between 
instruments  3(55  m)  and  7(15  m) .  Instrument  7  is  in  a  bottom  mixed 
layer  for  60  per  cent  of  the  record,  while  3  never  is.  The  coherence 
for  both  components  is  high  at  tidal  and  inertial  frequencies  and  drops 
to  below  the  95  per  cent  significance  level  above  the  local  Vaisala 
frequency.  The  only  significantly  non-zero  phases  in  figure  4  occur 
at  inertial  and  tidal  frequencies.  At  the  inertial  frequency  the 
boundary  layer  leads  the  interior  by  about  25°.  Similar  results  were 
found  by  Hayes  (1980)  and  Kundu  (1976) .  Chapter  7  shows  that  these 
phases  are  much  larger  than  would  be  expected  for  internal  waves  alone 
and  must  therefore  be  attributed  to  boundary  layer  effects. 

Significantly,  the  vertical  temperature  coherences  (not  shown)  in 
this  experiment  are  below  the  95  per  cent  significance  level  at  most 
frequencies,  for  separations  greater  than  30  m.  This  is  in  marked 
contrast  to  measurements  in  the  oceanic  interior,  where  the  temperature 
coherences  are  much  higher  than  velocity  coherences  (Briscoe,  1975) . 

This  is  clearly  due  to  the  strong  temperature  signals,  associated  with 
the  bottom  mixed  layer,  which  are  not  present  in  the  interior. 

Each  mean  spectrum  shown  in  figures  2  and  3  is  the  average  of 
measurements  from  a  particular  instrument  at  a  fixed  height  above  the 
bottom.  The  mixed  layer,  however,  varies  greatly  in  height,  so  that  any 
particular  near  bottom  instrument  is  sometimes  within  the  mixed  layer 
and  sometimes  above  it.  Averaging  over  measurements  at  a  fixed  height 
obscures  any  structure  associated  with  the  mixed  layer.  An  alternative 
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Figure  4 

Cross  spectra  between  east  (heavy  line)  and  north  (light  line)  velocity 
comi>onents  of  15  m  and  55  m  instruments,  computed  as  figure  1.  The  95% 
level  of  no  significance  and  95%  confidence  limits  for  zero  phase  are 
shown  as  dotted  lines  (Koopmans,  1974) .  The  15  m  level  leads  the  85  m 
level  at  the  inertial  frequency. 


approach,  which  will  be  employed  here,  is  to  filter  the  observed  veloci¬ 
ties  to  form  the  time  series  of  velocity  within  a  particular  spectral 
band.  The  variations  in  this  filtered  velocity  are  then  compared  with 
the  observed  density  structure,  so  that  features  of  the  velocity  field 
which  are  correlated  with  the  density  structure  may  be  identified.  This 
technique  has  been  applied  to  frequency  bands  spanning  the  superinertial 
range  of  frequencies  shown  in  figure  2.  Significant  structure  is  found 
only  for  high  frequency  (>  1  cph)  and  near  inertial  bands,  as  might  be 
expected  from  the  mean  spectra.  The  results  for  these  two  bands  are  dis¬ 
cussed  in  the  following  two  chapters. 


Chapter  4.  Observations  of  high  frequency  velocities 

The  benthic  boundary  layer  is  believed  to  be  turbulent.  An  in¬ 
creasing  level  of  high  frequency  velocity  fluctuations  is  seen  in  this 
experiment  near  the  bottom  (figure  2).  In  this  chapter  these  fluct¬ 
uations  are  shown  to  be  consistent  with  the  expected  energy  level  of 
boundary  layer  turbulence,  and  their  intensity  is  correlated  with  the 
observed  density  and  lower  frequency  velocity  fields. 

East  and  north  velocities  from  each  instrument  were  high  pass 
filtered  with  an  11  element  smooth  non-recursive  filter  constructed  by 
an  algorithm  described  in  Hamming  (1977) .  The  frequency  response  of 
this  filter  is  shown  in  figure  2.  Large  spikes  in  the  filtered  signal 
are  caused  by  instrument  rotor  stalls.  These  spikes  have  been  eliminated 
in  the  processing.  The  velocity  signals  processed  in  this  way  will  be 
called  high  pass  velocity. 

The  velocity  measured  by  the  VACM  sensors  (rotor  and  vane)  is 
vector  averaged  by  the  VACM  electronics  over  a  7.5  minute  interval.  The 
recorded  velocity  is  thus  a  low  pass  filtered  version  of  the  velocity 
measured  by  the  sensors.  The  combination  of  this  VACM  low  pass  filter 
and  the  digital  high  pass  filter  described  above  yields  a  band  pass 
filter  with  a  maximum  energy  transmission  of  .6  and  an  equivalent 
50  per  cent  energy  transmission  pass  band  from  2  to  5  cph.  Assuming  the 
VACM  sensors  to  be  perfect,  the  transfer  function  between  the  oceanic 
velocity  and  the  computed  high  pass  velocity  is  described  by  this  band 
pass  filter. 
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Observed  high  pass  velocities 


Figures  5,  6  and  7  display  the  high  pass  velocity  signal  for  3 
periods.  Various  features  of  the  observed  high  pass  velocity  are  illus¬ 
trated  in  these  figures  and  are  discussed  below.  The  complete  record  of 
high  pass  energy  is  displayed  in  figure  8c.  The  corresponding  potential 
temperature  and  1  day  mean  speed  records  are  shown  in  figures  8b  and  8a. 

A  thick  mixed  layer:  Figure  5a  shows  the  high  pass  north  velocity  for 
each  instrument.  The  amplitude  of  the  high  pass  fluctuations  is  clearly 
modulated  at  slower  time  scales.  Figure  5b  displays  contours  of  filtered 
high  pass  energy,  as  a  function  of  time  and  height  off  the  bottom.  The 
energy  time  series  are  calculated  from  the  high  pass  east  and  north 
velocity  time  series.  These  are  then  filtered  with  a  half  cosine  filter 
of  half  width  3  hours  to  form  the  high  pass  energy  time  series  displayed 
in  figure  5b.  The  correspondence  between  regions  of  large  high  pass 
energy  in  figure  5b  and  large  amplitude  fluctuations  in  figure  5a  should 
be  clear. 

Figure  5c  shows  2  m°C  isotherms  of  potential  temperature  as  a  fun¬ 
ction  of  time  and  height  off  the  bottom.  The  potential  temperature  here, 
and  in  all  subsequent  figures,  is  filtered  with  a  half  cosine  filter  of 
half  width  3  hours,  the  same  filter  used  on  the  high  pass  energy  in 
figure  5b.  A  45  m  thick,  well  defined  bottom  mixed  layer  is  present  from 
July  26  to  August  1.  The  lower  frequency  velocity  during  this  period 
is  approximately  the  same  at  all  levels  and  consists  or  a  steady  7  cm/s 
current  with  approximately  1  cm/s  internal  wave  fluctuations  (figure  8a) . 


19 


Figure  5 


High  pass  velocity  during  a  period  of  thick  mixed  layer 

a)  High  pass  north  velocity  for  all  7  instruments.  Frequency  response 
of  high  pass  filter  is  shown  in  figure  2.  Data  from  each  instru¬ 
ment  is  plotted  so  that  zero  velocity  corresponds  to  the  height 

of  that  instrument  off  the  bottom  as  shown  on  the  left  hand  axis. 
The  data  at  75  m  is  formed  by  averaging  the  65  m  and  85  m  data. 

b)  Contours  of  filtered  high  pass  energy,  plotted  as  a  function  of 
time  and  height  off  the  bottom  on  the  same  scale  as  a) .  The 
energy  is  formed  from  the  individual  high  pass  velocity  components 
and  then  filtered  with  a  half  cosine  filter  of  half  width  3  hours. 
Notice  how  clouds  of  high  pass  energy  intermittently  fill  the 
mixed  layer  shown  in  c) . 

c)  Contours  of  potential  temperature  filtered  as  in  b)  plotted  on 
same  scale. 

d)  Least  squares  ellipses  fit  to  successive  24  hour  pieces  of  high 
pass  velocity.  The  ellipses  are  rotated  so  that  high  pass 
velocity  fluctuations  parallel  to  the  24  hour  mean  velocity  are 
horizontal  on  the  figure;  velocity  fluctuations  perpendicular  to 
the  24  hour  mean  velocity  are  vertical.  The  size  and  orientation 
of  these  ellipses  show  the  energy  and  directionality  of  the  high 
pass  velocity,  respectively. 
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Comparing  the  high  pass  energy  seen  in  figure  5b  with  the  density  struc¬ 
ture  shown  in  figure  5c  reveals  clouds  of  high  pass  energy  extending  from 
the  bottom  and  intermittently  filling  the  mixed  layer.  This  looks  like 
bottom  generated  turbulence.  It  is  more  intense  near  the  bottom,  and  is 
present  only  when  the  mixed  layer  is  present.  Below,  the  energy  and 
frequency  of  these  fluctuations  is  shown  to  be  similar  to  that  expected 
for  bottom  turbulence. 

Strong  density  gradients:  Figures  6b  and  6c  show  high  pass  energy  and 
potential  temperature  contours  as  in  figures  5b,  c.  The  density  struc¬ 
ture  is  more  complex  here  with  a  frontlike  feature  on  June  4  and  5,  a 
shallow  mixed  layer  from  June  5  to  9  and  a  thick  mixed  layer  starting  on 
June  9.  Note  that  the  potential  temperature  contour  interval  is  5  m°C 
as  opposed  to  the  2  m°C  in  figure  5.  The  density  gradients  above  the 
mixed  Layer  are  much  stronger  here  than  in  figure  5.  The  mean  speed  is 
also  much  lower  during  this  period.  Figure  6a  marks  times  during  which 
the  VACM  rotor  did  not  move.  During  these  periods  the  calculated  high 
pass  velocity  is  zero. 

Figure  6b  shows  the  high  pass  energy  for  this  period.  Intermittent 
near  bottom  maxima  of  high  pass  energy,  similar  to  those  in  figure  5b, 
are  seen  here  within  the  thick  mixed  layer  starting  on  June  9.  Here, 
however,  they  are  much  less  energetic,  and  do  not  extend  to  the  mixed 
layer  top.  Unlike  in  figure  5b,  much  of  the  high  pass  energy  is  here  co¬ 
herent  with  tho  regions  of  strong  temperature  gradient;  the  front  on 
June  4,  5  and  the  mixed  layer  cap  starting  on  June  7  clearly  correspond 
with  regions  of  large  high  pass  energy.  The  high  pass  energy  associated 
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Figure  6 

High  pass  velocity  during  a  period  of  strong  density  gradients. 

a)  Regions  of  rotor  stall  for  each  instrument,  plotted  as  in  figure  5a 
Dark  regions  indicate  a  stall. 

l>)  As  figure  5b. 

c)  As  figure  5c,  but  5m°C  contour  interval. 


with  these  temperature  gradients  is  distinct  from  the  near  bottom,  pre¬ 
sumably  turbulent,  high  pass  energy.  It  may  result  from  advection  of 
small  scale,  low  frequency  shear,  associated  with  the  temperature 
gradients,  by  the  low  frequency  horizontal  current.  Significantly,  the 
interior  homogeneous  region  extending  from  June  5  to  8,  shows  little 
high  pass  energy.  If  this  is  a  detached  bottom  mixed  layer,  as  hypo¬ 
thesized  by  Armi  and  D'Asaro  (1980)  it  clearly  does  not  show  the  same 
level  of  high  pass  energy  as  the  bottom  mixed  layers.  This  suggests  that 
it  is  no  longer  actively  turbulent. 

A  complex  region:  Figure  7  sho-  s  the  high  pass  energy  and  potential 
temperature  structure  during  a  third  period,  characterized  by  rapid 
changes  in  the  mixed  layer  characteristics  and  a  large  (10  cm/s)  mean 
current.  During  this  period  interior  homogeneous  layers  merge  with  bot¬ 
tom  mixed  layers  on  3  occasions:  June  26,  June  29  and  July  2.  Mixed 
layer  detachment  may  be  occurring  here  as  hypothesized  by  Armi  and 
D'Asaro  (1980). 

The  high  pass  energy  structure  shows  many  features  seen  previously. 
The  bottom  intensified  maxima  on  June  27,  July  2,  and  July  3  appear  sim¬ 
ilar  to  those  in  figure  5;  the  maxima  on  July  3  extends  to  the  mixed 
layer  top  50  m  above  the  bottom.  The  region  of  strong  temperature  gra¬ 
dient  from  June  23  to  June  26  is  accompanied  by  regions  of  large  high 
pass  energy,  as  in  figure  6.  The  numerous  interior  homogeneous  layers 
show  a  generally  low  level  of  high  pass  energy.  The  maxima  in  high  pass 
energy  on  June  26,  however,  occurs  within  a  bottom  mixed  layer,  which  is 
continuous  with  an  interior  layer.  This  feature  is  unusual  in  that  it 
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Figure  8 


Mean  velocity,  density,  high  pass  velocity  and  near  inertial  velocity 
structure  for  entire  record.  All  figures  plotted  with  same  vertical 
and  horizontal  scales.  Differences  in  start  and  stop  times  are  due  to 
differences  in  averaging. 

a)  Left  axis:  1  day  mean  speed  at  15  m.  Right  axis:  Turbulent 
Ekman  layer  height,  .4u*/f,  computed  using  u*  =  S/30. 

S  is  1  day  mean  speed. 

b)  Potential  temperature  contours  as  in  figure  6c. 

c)  High  pass  energy  contours  as  in  figure  5b. 

d)  High  pass  velocity  ellipses  as  in  figure  5d. 

e)  Contours  of  near  inertial  clockwise  energy  computed  as  described 
in  the  text.  The  spectral  response  of  the  near  inertial  filter 
is  shown  in  figure  2.  A  nearly  independent  estimate  is  computed 
every  60  hours. 

f)  Contours  of  near  inertial  phase  of  each  instrument  relative  to 
the  phase  of  the  vertical  mean  near  inertial  velocity.  Clockwise 
leading  phases  are  shaded.  Phase  contours  are  not  shown  when 
near  inertial  energy  is  less  than  .25  cm2/s2. 
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is  not  bottom  intensif red ,  but  occurs  within  a  well  mixed  region.  A 
separated  turbulent  baundary  layer  would  have  such  a  mid -water  maximum 
in  turbulent  kinetic  energy. 

Expected  Turbulence  Properti es 

The  above  observations  reveal  two  major  kinds  of  high  pass 
velocity  fluctuations,  fluctuations  occurring  within  regions  of  strong 
temperature  gradient,  and  bottom  intensified  fluctuations  occurring 
within  the  bottom  mixed  layer.  This  second  variety  of  fluctuations  looks 
like  bottom  generated  turbulence.  In  this  section  the  expected  pro¬ 
perties  of  such  bottom  generated  turbulence  are  defined  and  a  comparison 
with  the  observed  properties  is  made. 

The  turbulence  in  the  benthic  boundary  layer  is  probably  generated 
by  a  mechanism  similar  to  that  operating  in  a  two  dimensional  laboratory 
turbulent  boundary  layer,  such  as  can  be  produced  in  a  wind  tunnel  or 
flume.  Such  boundary  layers  have  been  extensively  studied  and  their 
properties  are  well  known.  The  extrapolation  of  these  measurements  to 
the  benthic  boundary  layer  will  define  the  expected  properties  of  the 
boundary  layer  turbulence. 

Friction  velocity:  The  turbulent  velocity  fluctuations  within  the  boun¬ 
dary  layer  are  expected  to  scale  with  u*,  the  friction  velocity,  where 
pu*7  equals  the  wall  drag.  Since  u*  is  not  measured  in  this  experi¬ 
ment  it  must  be  estimated  from  the  measured  velocity  at  15  m.  The  impor¬ 
tant  parameters  in  this  estimate  are  the  bottom  roughness  and  the 
Reynold’s  number,  Re  =  1^6/v.  An  interior  velocity  =  10  cm/s  and 


a  boundary  layer  height  of  6  =  20  m  yield  Re  =  2  x  io6.  For  a  smooth 
wall  turbulent  boundary  layer  at  this  Re,  a  value  of  u^/u*  3  35  is 
predicted  using  the  Millikan  drag  law  (Millikan,  1039) .  Bottom  photo¬ 
graphs  of  the  Hatteras  Abyssal  Plain  show  a  flat  bottom  with  occasional 
1  cm  bumps  (Biscaye  and  Eittreim,  1977).  A  sand  grain  roughness  of 
i  =  1  cm,  which  certainly  overestimates  the  effect  of  these  bumps, 
yields  a  roughness  Reynold's  number  u*£/v  of  approximately  30,  which 
is  transitional  between  smooth  and  rough  flow  (Hinze,  1959,  p.  486)  and 
predicts  U,/u*  y  29.  Csanady's  (1967)  Ekman  drag  formula  with 

z  =  v/u*  yields  U  /u*  =34.  A  value  of  U,/u*  =  32  will  be  used 
o  1  1 

below. 

Frequency  and  energy:  Figure  9a  shows  the  measured  turbulent  velocity 
fluctuation  intensity  through  a  laboratory  turbulent  boundary  layer 
(Klebaooff,  1954).  The  turbulent  velocity  fluctuations  scale  with  u* 
and  become  small  above  ^gg»  the  boundary  layer  thickness.  In  the  outer 
part  of  the  boundary  layer  the  horizontal  kinetic  energy  due  to  these 
turbulent  velocity  fluctuations  is  approximately  u*2.  Using 
Uj/u*  =  32,  this  yields  u*2  =  10  3  U^2. 

The  high  pass  energy  is  sensitive  only  to  velocity  fluctuations 
within  a  particular  range  of  frequencies.  Do  the  expected  turbulent 
velocit,  fluctuations  have  energy  in  this  frequency  range?  Figure  9b 
shows  the  spectra  of  downstream  and  crossstream  velocity  fluctuations 
in  the  outer  part  of  a  laboratory  turbulent  boundary  layer  (Bradshaw, 
1966) .  Velocity  spectra  in  the  outer  part  of  a  laboratory  turbulent 
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Ekman  layer  are  similar  (Caldwell  et  al.  1972) .  Plotted  in  the  same 
non-dimensional  coordinates  is  the  high  pass  energy  central  frequency, 

3.5  cph,  for  a  range  of  mean  speeds,  5-10  cm/s,  and  a  range  of  6 
10-20  m.  The  high  pass  energy  central  frequency  is  seen  to  be  near  the 
spectral  peak  of  the  expected  turbulent  velocity  fluctuations.  The  high 
pass  energy  should  thus  be  an  excellent  indicator  of  the  presence  of 
turbulent  boundary  layer  velocity  fluctuations. 

The  half  power  bandwidth  of  the  high  pass  velocity  transfer  func¬ 
tion  is  also  shown  in  figure  9b.  Roughly  20  per  cent  of  the  turbulent 
horizontal  kinetic  energy  shown  in  figure  9b  is  contained  within  this 
bandwidth.  Since  this  is  a  50  per  cent  transmission  bandwidth  a 
measurement  of  high  pass  energy  in  the  outer  part  of  a  turbulent  boun¬ 
dary  ls  expected  to  yield  approximately  10  per  cent  of  the  total  tur¬ 
bulent  horizontal  kinetic  energy  or  (.01  U^)2. 

Stratification:  The  mixed  layers  observed  in  this  experiment  commonly 

show  a  slightly  stable  density  stratification,  measured  by  potential  tem¬ 
perature.  Stable  density  stratif ication  can  strongly  influence  a  tur- 
bui  1*11*  txiundary  layer  and  if  sufficiently  strong  can  suppress  the  tur¬ 
bulence  completely  (Turner,  1973;  Arya,  1972).  Assuming  the  energetic 
eddies  in  the  outer  part  of  the  boundary  layer  to  be  of  the  same  scale 
as  the  boundary  layer,  II  =  10  m  and  that  they  have  an  energy  density 
u*'  ,  Pi  7  cm,  .  mves  a  kinetic  energy  of  u*  H  y  48  ergs/cm  . 

Assuming  that  these  eddies  will  completely  mix  any  existing  stratifica- 

Ap  H2  „  . 

Mon  on  *  he  scale  II,  this  will  require  energy  g^  jj.  Using 

;  —  . i  ••  and  H  =  10  m,  4  ergs/cm2  are  required  to  mix  a  tem- 


perature  gradient  of  .5  m°C/10  m,  typical  of  a  mixed  layer,  while 
80  ergs/cm2  are  required  to  mix  the  typical  10  m°C/10  m  temperature 
gradient  above  the  mixed  layer.  Thus  the  typical  stratification  within 
the  mixed  layer  can  be  mixed  by  the  boundary  layer  turbulence  without 
suppressing  this  turbulence.  Boundary  layer  turbulence  in  the  strati¬ 
fied  fluid  above  the  mixed  layer,  however,  will  be  strongly  damped  by 
the  stratification.  Velocity  fluctuations  in  the  stratified  region  are 
thus  likely  to  be  wavelike,  rather  than  turbulent  (Csanady,  1978; 

Turner,  1973;  Piat  and  Hopfinger,  1980).  Boundary  layer  turbulence  is 
therefore  expected  only  within  the  bottom  mixed  layer. 

Is  it  turbulence? 

k 

Figure  10  shows  scatter  plots  of  (high  pass  energy) "  versus  speed 
at  85  m  (10a)  and  15  m  (10b)  for  the  entire  92  day  record.  Speed  is 
here  formed  from  the  east  and  north  velocity  components  filtered  with  a 
half  cosine  filter  of  3  hour  half  width.  The  15  m  instrument  is  commonly 
in  a  mixed  layer,  while  the  85  m  instrument  never  is.  At  85  m  the  high 
pass  energy  displays  a  constant  level  of  approximately  (.04  cm/s)2 inde¬ 
pendent  of  speed  except  at  low  speeds  where  the  finite  threshold  of  the 
VACM  is  im(>ortunt.  At  15  m  the  fluctuation  intensity  is  much  higher, 
approximately  .1  em/s  at  a  speed  of  10  cm/s  and  is  roughly  correlated 
with  speed.  At  15  m  the  high  pass  energy  level  is  consistent  with  the 
expected  properties  of  boundary  layer  turbulence,  i.e.,  the  velocity 
fluctuations  scale  with  speed,  and  have  a  magnitude  of  .01U.  This  is 


this  instrument  is  commonly  within  mixed  l-'.yeis. 


The  observed  high  pass  velocity  fluctuations  appear  to  be  due 
primarily  to  boundary  layer  turbulence.  This  is  supported  by  the  fol¬ 
lowing  observations:  The  high  pass  velocity  is  sensitive  to  velocity 
fluctuations  near  the  expected  spectral  peak  of  boundary  layer  tur¬ 
bulence  (figure  9b).  The  observed  high  pass  energy  level  is  roughly  cor¬ 
related  with  the  lower  frequency  velocity  (figure  10b)  as  expected  for 
boundary  layer  tv-bulence;  the  high  pass  energy  level  is  consistent  with 
the  expected  turbulent  energy  level.  Boundary  layer  turbulence  should 
be  confined  to  the  mixed  layer,  and  be  more  energetic  near  the  ocean  bot¬ 
tom.  The  high  pass  energy  in  figure  5  is  certainly  of  this  character, 
as  is  much  of  the  high  pass  energy  in  figures  6  and  7. 

Sources  of  high  frequency  velocity  fluctuations  other  than  boun¬ 
dary  layer  turbulence  are  also  apparent  in  the  data.  The  velocity  fluc¬ 
tuations  associated  with  the  large  temperature  gradients  in  figures  6 
and  7  are  one  example.  The  background  level  of  high  pass  energy  in 
fiqure  10a  suggests  an  instrumental  noise  source.  Simulations  of  VACM 
round  off  errors  (McCullough,  1975)  yield  an  intrinsic  high  pass  velo¬ 
city  noise  level  of  about  .03  cm/s.  This  is  comparable  to  the  back¬ 
ground  energy  level  in  figure  10a.  Other  instrumental  sources  of  high 
pass  energy,  from  mooring  motion  or  imperfect  sensor  response,  are 
possible.  These  sources,  however,  should  be  sensitive  only  to  the  velo¬ 
city  field,  not  to  the  boundary  layer  density  structure.  The  observed 
fluctuations  show  a  clear  correlation  with  the  boundary  layer  density 
structure.  The  high  pass  velocity,  especially  that  part  occurring  with- 
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in  bottom  mixed  layers,  will  thus  be  interpreted  as  due  primarily  to 
boundary  layer  turbulence. 

The  high  pass  velocity  is  a  measure  of  boundary  layer  turbulence. 
Unfortunately,  it  is  a  quite  limited  one.  It  is  sensitive  only  to  velo¬ 
city  fluctuations  in  a  narrow  frequency  range,  which  is  estimated  to  con¬ 
tain  only  10%  of  the  total  turbulent  kinetic  energy.  Changes  in  either 
the  frequency  or  energy  of  the  turbulent  velocity  fluctuations  will  re¬ 
sult  in  a  change  in  the  high  pass  energy,  so  that  changes  in  high  pass 
energy  do  not  necessarily  correspond  to  changes  in  the  turbulent  energy. 
The  high  pass  velocity  is  also  a  noisy  measure  of  boundary  layer  tur¬ 
bulence  with  a  significant  high  pass  signal  coming  from  other  oceanic 
velocity  fluctuations  (figure  6)  and  instrumental  noise.  The  high  pass 
velocity  can  therefore.'  not  be  used  to  detail  the  structure  of  boundary 
layer  turbulence,  but  only  to  indicate  its  approximate  strength  as  a 
function  of  time.  The  measurements  are  further  limited  by  the  VACM 
locations  to  between  15  and  85  m  aoove  the  bottom.  Weatherly’s  (1972) 
observation:.;  of  the  boundary  layer  with  a  similar  mean  speed  indicate 
that  thus  is  significantly  above  the  logarithmic  region  of  the  boun¬ 
dary  layer.  Kememboring  all  of  the  above  limitations,  the  measurements 
of  high  pass  energy  provide  a  useful  measure  of  the  turbulence  struc¬ 
ture  across;  the  mixed  layer,  which  can  be  compared  with  other  measure¬ 
ments  and  theoretical  constructs. 


Observations  of  boundary  layer  turbulence 


The  high  pass  velocity  shows  several  interesting  aspects  of  the 

boundary  layer  turbulence.  The  turbulence  appears  to  intermittently 

fill  the  entire  bottom  mixed  layer.  This  is  seen  clearly  in  figure  5, 

starting  on  July  26;  in  figure  7  on  July  3,  and  in  figure  8  on  May  18-21, 

May  30,  July  11,  and  August  13-18.  In  figure  8a  the  turbulent  Ekman 

layer  height  6  =  .4u*/f  is  plotted  with  u*  =  S/30,  S  being  the 

EK 

speed  computed  from  the  1  day  mean  velocity  at  15  m.  6  is  generally 

EK 

less  than  15  m  and  usually  is  less  than  the  mixed  layer  height  (Armi  and 

Millard,  1976).  For  example  in  figure  5  6  „  y  15  m,  while  the  mixed 

EK.  n, 

layer  height  and  the  penetration  height  of  the  boundary  layer  turbulence 
is  approximately  40  m.  Other  examples  can  be  seen  in  figure  7  on  July  3 
and  in  figure  8  on  May  18-21  and  August  13-18,  while  a  counter  example 
can  be  seen  in  figure  6  starting  on  June  9.  These  observations  indicate 
that  the  boundary  layer  turbulence  commonly  extends  above  a  turbulent 
Ekman  layer  height,  in  contrast  to  the  expected  properties  of  a  turbulent 
Ekman  layer.  These  observations  are  consistent  with  previous  obser¬ 
vations  of  Radon-222  (Sarmiento,  1978)  in  this  area,  which  indicate  that 
turbulent  mixing,  on  the  time  scale  of  a  few  days,  penetrates  to  the  top 
of  the  mixed  layer,  and  above  a  turbulent  Ekman  layer  height. 

Fiqures  5d  and  8d  display  the  directionality  of  the  high  pass  hori¬ 
zontal  velocity.  For  each  24  hour  period  the  best  least  squares  ellipse 
is  fit  to  the  two  components  of  high  pass  velocity  (Caiman,  1978)  .  These 
ellipses  are  plotted  for  each  day  and  rotated  so  that  a  vertical  line 
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corresponds  to  high  pass  velocity  fluctuations  oriented  perpendicular 
to  the  14  hour  mean  velocity,  while  a  horizontal  line  corresponds  to 
fluctuations  oriented  parallel  to  the  24  hour  mean  velocity.  Figure  5d 
clearly  shows  a  concentration  of  high  pass  energy  within  the  mixed 
layer  and  near  the  bottom.  The  velocity  fluctuations  are  seen  to  align 
themselves  roughly  45°  clockwise  of  the  mean  velocity  whenever  the 
bottommost  instrument  is  within  a  mixed  layer.  A  preference  for  this 
general  direction  can  also  be  seen  in  figure  8d  (May  18-21,  May  28, 

July  5-10,  August  12-16).  Assuming  that  these  turbulent  fluctuations 
are  generated  by  the  lower  frequency  mean  velocity,  this  asymmetry 
between  the  direction  of  the  mean  velocity  and  the  orientation  of  the 
fluctuations,  suggests  that  rotation,  which  enters  the  equations  of 
motion  in  an  asymmetrical  way,  must  be  important  in  the  dynamics  of  the 
boundary  layer. 

The  discussion  of  these  observations  is  postponed  to  chapter  6, 
so  that  the  observations  of  the  near  inertial  motions,  which  yield  ad¬ 
ditional  information  on  the  boundary  layer  turbulence,  may  be  described. 


Chapter  5.  Observations  of  near  inertial  velocities 

The  average  horizontal  kinetic  energy  spectra  shown  in  figure  2 
indicate  less  near  inertial  energy  in  the  bottommost  instruments  than 
in  the  upper  instruments.  Similar  near  bottom  decreases  in  near  inertial 
energy  have  been  observed  by  Kundu  (1976) ,  Hayes  (1980)  and  Leaman  (1976) 
and  can  clearly  be  seen  in  figure  5  of  Weatherly  and  Wimbush  (1980) .  In 
the  following  section  the  relation  of  this  decrease  to  the  bottom  mixed 
layer  structure  is  examined. 

The  92  day  velocity  record  from  each  instrument  was  broken  into  36, 
120  hour  long  50%  overlapping  pieces.  The  near  inertial  clockwise  and 
anticlockwise  components  of  velocity  in  each  piece  were  found  using  a 
demodulation  filter  with  center  frequency  of  1/24  cph  and  a  4-term 
Blackman-Harris  window  described  by  Harris  (1978) .  The  spectral  re¬ 
sponse  of  this  window  is  shown  in  figure  2.  Near  the  center  frequency 
it  has  a  shape  similar  to  that  of  a  Hanning  window  with  a  72  hour  piece. 
The  side  lobes,  however,  are  much  lower.  Several  other  windows  have 
been  used  and  none  of  the  results  shown  below  depend  critically  on  the 
particular  window  used. 

Clockwise  velocities 

Most  of  the  near  inertial  energy  is  in  clockwise  motions  (figure  3). 
The  analysis  of  the  near  inertial  velocity  will  thus  concentrate  on  the 
clockwise  component  of  motion.  Figure  8e  shows  contours  of  clockwise 
Lnertial  energy  as  a  function  of  time  and  height  off  the  bottom.  The 
near  inertial  energy  level  is  vertically  coherent,  but  varies  signifi- 
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cantly  with  time.  There  is  clearly  less  near  inertial  energy  near  the 
ocean  bottom,  as  is  seen  in  the  average  spectra  (figure  2) .  The  height 
at  which  the  energy  decreases  most  rapidly,  however,  varies  with  time. 
Comparing  8e  with  the  5  m°C  isotherms  of  potential  temperature  displayed 
in  figure  8b,  the  height  at  which  the  inertial  energy  decreases  most 
rapidly  is  seen  to  correlate  with  the  height  of  the  bottom  mixed  layer. 
For  example,  both  the  mixed  layer  height  and  the  height  of  maximum 
gradient  are  large  on  May  20  and  July  27-31.  On  June  21-29  and  August 
12-16  they  both  are  small.  The  mixed  layer  top  and  the  height  of  max¬ 
imum  gradient  both  decrease  20  m  between  July  6  and  July  13.  The  major 
exception  to  this  pattern  occurs  on  May  28,  when  a  maximum  in  inertial 
enerqy  is  seen  at  the  mixed  layer  top.  At  this  time,  however,  the  mean 
velocity  is  small  (figure  8a)  and  any  turbulent  boundary  layer  effects 
will  thus  be  small.  In  general  the  mixed  layer  exhibits  less  near 
inertial  energy  than  the  interior,  with  the  maximum  gradient  in  energy 
occurring  slightly  above  the  top  of  the  mixed  layer. 

Fiyure  8f  shows  contours  of  the  phase  of  the  clockwise  inertial 
velocity  at  each  level  relative  to  the  phase  of  the  mean  clockwise 
inertial  velocity.  Using  the  vertical  mean  as  a  reference  rather  than 
a  particular  level  allows  all  levels  to  be  treated  identically.  Regions 
of  negative  phase  lead  the  mean  and  are  shaded.  The  computed  phases 
are  not  shown  when  the  near  inertial  energy  is  less  than  .25  cm2/s2; 
they  are  quite  noisy  at  these  times. 
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The  computed  phases  clearly  show  the  15  m  instrument  leading,  on 
the  average,  relative  to  the  mean,  and  leading  almost  always  relative 
to  the  25  m  instrument.  The  typical  phase  difference  of  25°  is  sig¬ 
nificantly  larger  than  any  expected  instrumental  errors.  Unlike  the 
energy,  the  near  inertial  phase  exhibits  large  gradients  within  the  bot¬ 
tom  mixed  layer,  for  example  on  May  20-26,  June  3,  July  3  and  August  10- 
15.  Similar  gradients  do  not  generally  appear  above  the  mixed  layer. 

The  phases  computed  here  clearly  show  the  apparent  instrumental 
direction  errors  discussed  in  chapter  2.  Instruments  2  and  3  show  a 
consistently  positive  (anticlockwise)  phase,  while  instruments  1  and  4 
show  a  consistently  negative  (clockwise)  phase.  This  effect  is  espec¬ 
ially  striking  in  late  June  and  early  July  in  the  long  tongues  of 
negative  phase  at  45  m  (instrument  4) . 

The  phase  structure  of  the  near  inertial  motions  shows  far  less 
correlation  with  the  density  structure  than  does  the  energy.  Although 
this  may  reflect  the  physics  of  the  boundary  layer  it  may  also  be  due 
to  the  much  higher  noise  level  in  the  phase  measurement.  Spectral 
leakage  from  adjacent  frequency  bands  affects  phase  more  strongly 
than  energy  since  these  other  frequencies,  by  definition,  have  rapidly 
varying  phase  with  respect  to  the  inertial  frequency.  The  spectral 
leakage  is  also  more  severe  for  phase,  being  proportional  to  the  filter 
transmission,  rather  than  the  transmission  squared  as  for  energy. 

Finally  the  phase,  unlike  horizontal  kinetic  energy,  is  sensitive  to 
instrumental  direction  errors. 
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Anticlockwise  velocities 

The  near  inertial  anticlockwise  motions  are  much  less  energetic 
than  the  clockwise  motions.  They  do  not  show  the  strong  correlation  with 
density  structure  found  for  the  clockwise  motions;  energy  and  phase 
digrams  similar  to  figures  8e,  f  are  therefore  not  shown  here. 

Simple  WKB  internal  wave  theory  predicts  a  ratio  T  of  clockwise 
to  anticlockwise  energy  of  T  =  (0  +  f)2/(0  -  f)2  for  0  f  (Fofonoff, 
1969).  Fu  (1980),  accounting  for  the  variability  of  f  with  latitude, 
finds  this  prediction  should  be  valid  for  a  >  l.Olf.  The  value  of  T 
observed  here  is  not  consistent  with  these  predictions.  For  example, 
figure  2  shows,  T  =  50,  at  a  =  f,  corresponding  to  a  WKB  prediction 
of  a  =  1.33f.  There  is  more  anticlockwise  energy  than  is  predicted. 

This  suggests  that  much  of  the  observed  anticlockwise  motion  is  not  due 
to  internal  waves.  Fu  (1980)  finds  a  similar  excess  of  near  inertial 
anticlockwise  energy  in  a  large  number  of  North  Atlantic  current  meter 
measurements . 

Anisotropy 

The  directionality  of  the  near  inertial  motions  can  be  examined 
by  computing  the  average  variance  ellipse  Following  Gonella  (1972) 
and  Caiman  (1978)  the  amplitude  and  phase  of 

C(0)  C(-a) 

K.  —  1 

|fi(a)  I  I c (-o) 

defines  the  stability  and  direction  of  the  average  ellipse.  C(O)  and 
C(-o)  are  the  clockwise  and  anticlockwise  amplitudes  for  the  near 
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inertial  band.  The  ellipse  stability  has  the  same  form  and  probability 
distribution  as  a  coherence;  if  it  is  significantly  different  from  zero, 
the  ellipse  is  significantly  anisotropic. 

Table  1  lists  the  ellipse  properties  computed  for  the  top  instru¬ 
ment  using  the  same  Fourier  transforming  procedure  used  to  make 
figure  8e,  f.  These  ellipse  parameters  are  computed  in  both  an  east/ 
north  (geographical)  coordinate  system,  and  with  each  Fourier  transformed 
piece  rotated  so  that  the  central  one  day  mean  velocities  in  each  piece 
are  parallel  (velocity  coordinate  system) .  The  one  day  mean  velocity 
points  roughly  east  for  most  of  the  record.  The  near  inertial  velocities 
show  a  small  but  significant  anisotropy,  in  either  coordinate  system, 
with  the  major  axis  of  the  ellipse  pointing  in  the  direction  of  the  mean 
velocity.  The  anisotropy  is  somewhat  stronger  in  the  velocity  coor¬ 
dinate  system. 

If  the  observed  velocity  fluctuations  are  interpreted  as  internal 
waves,  the  direction  of  the  velocity  ellipse  gives  the  direction  of 
wave  propagation  (Caiman,  1978) .  Here,  this  would  indicate  waves  pro¬ 
pagating  primarily  parallel  to  1  day  mean  velocity  vector.  This  inter¬ 
pretation  assumes,  however,  that  the  observed  anticlockwise  motions 
are  due  to  inter:, al  waves,  which,  as  discussed  above,  may  not  be  true. 


43 


Table  1 


Major  Axis  (cm/s) 
Minor  Axis  (cm/s) 

Direction  of  Major 

Ellipse  Stability 


Average  Near  Inertial  Ellipse 


Geographical 
Coordina tes 


.93 

.75 


Axis 


7°  from  East 


.36 


Velocity 

Coordinates 

.95 

.73 

3°  from  downstream 

.43 


95%  level  of  no 
significance  @  72  EDOF 


.29 


.29 


..ujnma 

Mutun.i  in  the  inertial -gravity  wave  frequency  band  show  a  stronq 
ouundary  Layer  signal  only  near  the  inertial  frequency.  The  near  inertial 
.  l.H.kwise  behavior  ls  summarized  in  figure  11.  Each  arrow  represents  the 
implitude  and  phase  of  the  near  inertial  clockwise  velocity.  The  energy 
i"  reuses  near  t  Lie  top  of  the  mixed  layer  with  little  phase  change. 

Within  the  mixed  Layer  the  phase  begins  to  lead,  with  the  phase  angle 
increasing  downward.  Extrapolating  to  the  bottom,  the  energy  should 
irop  to  zero,  with  an  unknown  phase.  These  variations  may  occur  sig¬ 
nificantly  above  a  turbulent  Ekman  layer  height. 

The  observed  near  inertial  anticlockwise  velocities  do  not  show 
a  similarly  simple  pattern  and  contain  more  energy  ttian  is  predicted  by 
simple  internal  wave  theory.  The  near  inertial  motions  are  slightly 
anisotropic  with  the  major  axis  of  the  mean  variance  ellipse  oriented 
parallel  with  the  1  day  mean  velocity. 
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The  one  dimensional  steady  boundary  layer  models  used  to  model  the 
benthic  boundary  layer  should  become  identical  to  the  turbulent  Ekman 
layer  model  in  the  limit  of  no  stratification.  In  particular  the  pre¬ 
dicted  thickness  should  become  the  turbulent  Ekman  layer  height  in  this 
limit.  The  presence  of  stratification  should  decrease  the  predicted 
boundary  layer  height  (Businger  and  Arya,  1974).  The  observations  here, 
therefore,  show  a  boundary  layer  much  thicker  than  would  be  predicted 
by  steady,  o  dimensional  models. 

A  striking  feature  of  the  high  pass  energy  shown  in  figure  5  is 
the  highly  intermittent  character  of  the  boundary  layer  turbulence.  At 
times  the  turbulence  extends  throughout  the  mixed  layer,  while  a  few 
hours  later  it  is  completely  absent.  These  variations  occur  on  tidal 
and  inertial  time  scales.  This  suggests  that  the  observed  fluctuations 
in  the  boundary  layer  may  be  driven  by  the  oceanic  internal  wave  field, 
which  has  similar  time  scales. 

A  laminar  Ekman  layer,  when  subject  to  a  clockwise  rotating 
interior  flow  of  frequency  0,  increases  in  height  by  the  factor 
| l  -  0/f I  leading  to  an  infinite  height  for  inertial  forcing 
(Greenspan,  1968) .  Figure  3  shows  that  near  inertial  clockwise 
motions  are  present  in  a  significant  amplitude.  This  suggests  that 
they  may  thicken  the  boundary  layer.  Weatherly  et^  aJU  (1980)  using  a 
turbulence  closure  model  to  simulate  the  boundary  layer  in  this  region 
find  that  the  addition  of  inertial  forcing  results  in  a  doubling  of 
the  boundary  layer  height.  This  mechanism  may  explain  the  thickness  of 
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;he  observed  boundary  layei .  It  does  not  address  the  question  of  inter- 
mittency . 

Turbulent  boundary  layers  are  known  to  be  quite  sensitive  to  the 
acceleration  or  deceleration  of  their  driving  flow.  For  two-dimensional 
boundary  layers  deceleration  results  in  a  decreased  drag  coefficient,  a 
more  rapid  growth  rate,  higher  turbulence  levels,  and  if  continued  suf¬ 
ficiently  long,  separation  (Yaglom,  1979;  Goldstein,  1965) .  It  seems 
plausible,  therefore,  that  the  observed  intermittency  in  the  benthic 
boundary  layer,  is  driven  by  the  intermittent  internal  wave  deceleration 
of  the  boundary  layer.  Application  of  a  simple  scaling  theory  (D’Asaro, 
1980)  to  the  benthic  boundary  layer  predicts  that  the  observed  internal 
wave  decelerations  should  lead  to  a  significant  increase  in  boundary 
layer  thickness  on  a  few  hours  time  scale.  The  data,  however,  indicate 
no  significant  correlation  between  high  pass  energy  at  15  m  and  the  de¬ 
celeration  of  the  15  m  velocity  in  either  a  downstream  or  crossstream 
direction.  A  correlation  between  the  internal  wave  vertical  velocity 
and  Ixiundary  layer  thickness  might  also  be  expected.  Using  the  35  m 
temperature  variations  to  compute  the  vertical  velocity,  no  correlation 
is  found.  Thus  although  the  observed  boundary  layer  fluctuations  may 
be  driven  by  the  internal  wave  field,  no  direct  evidence  of  this  is 
found . 

The  possibility  remains  that  the  observed  boundary  layer  height 
and  intermittency  are  not  due  to  the  unsteadiness  of  the  interior 
velocity,  but  is  instead  an  intrinsic  feature  of  the  benthic  boundary 
layer  driven  by  a  steady  velocity.  "Bursting"  intermittency,  with  a 
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period  of  roughly  2.5  6/u  is  commonly  observed  in  the  outer  part  of  tur 
bulent  boundary  layers  (Kovasznay  et  al.,  1970;  Falco,  1977)  where  - 
is  the  thickness  of  the  boundary  layer.  Using  6  =  20  m,  IJ  =  7  cm/s, 
a  "bursting"  period  of  less  than  an  hour  is  expected.  This  is  much 
shorter  than  the  observed  intermittency  period.  The  observed  intermit- 
tency  is  thus  probably  not  due  to  "bursting". 

A  laminar  Ekman  layer  is  unstable  at  sufficiently  high  Reynold's 
number.  Could  the  boundary  layer  fluctuations  observed  here  be  due  to 
a  similar  type  of  secondary  circulation  within  the  turbulent  Ekman  layer 
as  postulated  by  Brown  (1970)?  The  calculated  instabilities  for  a  lami¬ 
nar  Ekman  layer  (Lilly,  1966;  Faller  and  Kaylor,  1966)  have  a  frequency 
of  magnitude  much  larger  than  f(e.g.,  Caldwell  and  van  Atta,  1970, 
fig.  16)  which  changes  sign  very  near  the  fastest  growing  instability. 
Thus  the  predicted  frequency  of  the  fastest  growing  wave,  is  highly 
variable  depending  on  exactly  which  wave  is  chosen.  Most  likely,  how¬ 
ever,  the  frequency  is  much  larger  than  f,  and  thus  cannot  account 
for  the  boundary  layer  variations  observed  here.  It  should  be  noted 
that  Caldwell  and  van  Atta  (1970)  and  Deardorf f (1970)  see  no  evidence 
of  laminar  Ekman  layer  instabilities  at  the  high  Reynold's  numbers 
•ha r  ic ter i st  i c  of  the  turbulent  Ekman  layer. 

In  summary  then,  no  satisfactory  explantion  for  the  origin  of 
the  clouds  of  turbulence  seen  in  figure  5  has  been  found.  The  obser¬ 
vations  are  not  consistent  with  a  steady,  turbulent  Ekman  layer,  as  far 
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as  its  properties  are  known.  Although  internal  wave  forcing  may  cause 
the  boundary  layer  to  become  thicker  than  a  turbulent  Ekman  layer  height, 
no  correlation  between  the  observed  turbulence  fluctuations  and  observed 
internal  wave  fluctuations  is  found. 

The  observed  high  pass  and  near  inertial  velocities  both  show  a 
distinct  mixed  layer  signal  which  is  here  attributed  to  boundary  layer 
turbulence  within  the  mixed  layer.  Careful  examination  of  the  calculated 
high  pass  and  inertial  velocities  reveals,  however,  that  the  mixed 
layer  signal  for  both  extends  some  distance  into  the  stratified  fluid. 
This  suggests  that  the  boundary  layer  velocity  fluctuations  and  stresses 
extend  into  the  stratification,  perhaps  partially  as  evanescent  fluc¬ 
tuations  (Phillips,  1955)  or  high  frequency,  stress  carrying,  internal 
waves  (Csanady,  1978).  This  may  be  important  in  the  entrainment  of 
fluid  into  the  mixed  layer.  The  boundary  layer  model  of  Weatherly  and 
Martin  (1973)  exhibits  a  similar  behavior,  with  the  computed  turbulence 
intensity  and  the  diffusion  coefficients  diminishing  to  zero  only  above 
the  mixed  layer  top. 


51 


Chapter  7.  Internal  wave  interaction  with  the  boundary  layer 

During  the  last  decade  great  progress  has  been  made  in  describing 
the  space-time  structure  of  the  internal  wave  field  (Garrett  and  Munk, 
1979;  Gregg  and  Briscoe,  1979).  Still  elusive,  however,  is  a  clear 
understanding  of  the  dynamics  of  the  internal  wave  field  and  the  magni¬ 
tude  and  character  of  its  sources  and  sinks.  One  possible  location  of 
internal  wave  generation  and/or  absorption  is  the  benthic  boundary 
layer.  The  boundary  layer  is  known  to  have  characteristics  distinct 
from  the  oceanic  interior  (Bowden,  1978) ,  and  is  certainly  the  site  of 
considerable  energy  dissipation.  How  much  of  this  dissipated  energy 
comes  from  the  internal  wave  field?  Leaman  (1976)  estimates  a  net  down- 
ward  energy  flux  in  the  internal  wave  field  of  .2  -  .3  ergs/cm  /s  for 
a  4  day  period  in  the  Sargasso  Sea.  Muller  et  al.  (1978) ,  analyzing 
the  IWEX  experiment,  find  a  similar  asymmetry  confined  to  the  near  iner¬ 
tial  frequency  band.  Kundu  (1976)  finds  a  somewhat  larger  net  downward 
flux  in  the  near  inertial  band  on  the  Oregon  shelf.  One  interpretation 
is  that  this  net  downward  energy  flux  results  from  absorption  of  near 
inertial  waves  by  the  benthic  boundary  layer.  Fu  (1980),  however,  pro¬ 
poses  an  alternative  explanation  based  on  the  kinematic  properties  of 
surface  generated  near  inertial  waves.  Below  a  direct  estimate  of  the 
absorption  of  near  inertial  waves  by  the  benthic  boundary  layer  on  a 
flat  abyssal  plain  in  the  western  North  Atlantic  is  made. 
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Theoretical  approach 

The  near  inertial  observations  in  chapter  5  indicate  a  signifi¬ 
cant  interaction  between  the  near  inertial  motions  and  the  boundary 
layer.  In  the  analysis  presented  below,  these  motions  will  be  analyzed 
as  linear  internal  waves.  The  long  non-linear  interaction  times  com¬ 
puted  for  the  energy  containing  near  inertial  waves  (McComas,  1977) 
justifies  this  approach.  McComas'  calculations  are  based  on  thermo- 
cline  wave  scales.  In  the  deep  ocean  the  wave  scales  are  significantly 
larger,  and  the  non-linear  interactions  much  weaker. 

Consider  a  set  of  linear  internal  waves  incident  on  a  flat  bot¬ 
tomed  ocean,  with  a  thin  turbulent  bottom  boundary  layer.  These  waves 
are  superimposed  on  currents  with  much  slower  time  scales  here  called 
"mean"  currents.  In  the  absence  of  any  waves  the  mean  shear  within  the 
boundary  layer  is  balanced  by  the  turbulent  stress  divergence.  The 
energy  to  maintain  the  stresses  is  drawn  from  the  kinetic  energy  of  th 
mean  flow.  Any  wave  superimposed  on  this  balance  will  modulate  the 
stresses.  The  energy  to  maintain  the  stresses  may  now  come  partially 
from  the  wave  as  well  as  from  the  mean  flow.  If  there  are  many  waves, 
the  energy  may  be  drawn  from  some  unknown  combination  of  these  waves 
and  the  mean  flow.  Since  the  boundary  layer  is  highly  non-linear  one 
should  not  expect  the  energy  extracted  from  any  given  wave  to  depend 
solely  on  the  characteristics  of  that  wave;  it  may  also  depend  on  the 
other  waves  and  on  the  mean  flow. 
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The  following  analysis  is  designed  to  compute  the  absorption  of 
energy  from  linear  internal  waves  by  a  turbulent  boundary  layer.  The 
analysis  assumes  the  waves  to  be  large  compared  to  the  boundary  layer 
and  uses  the  observed  wave  motions  at  a  particular  frequency  to  compute 
the  turbulent  stresses  acting  on  these  motions.  No  assumption  concern¬ 
ing  the  source  of  these  stresses  is  made;  they  may  be  due  to  motions 
at  other  frequencies. 

Below,  the  relevant  internal  wave  and  boundary  layer  scales  will 
be  computed,  and  the  equations  of  motion  simplified  for  these  scales. 

The  resulting  equations  will  be  used  to  interpret  the  structure  shown  in 
figure  11;  the  first  qualitatively,  and  then  quantitatively  by  computii  g 
the  energy  absorption  from  the  near  inertial  waves. 

Wave  scales 

The  vertical  scales  of  internal  waves  in  the  thermocline  away  from 
the  inertial  frequency  are  well  known  (Garrett  and  Munk,  1979) .  The 
vertical  scales  of  near  inertial  waves,  particularly  in  the  deep  ocean, 
are  less  well  characterized.  Munk  and  Phillips  (1968)  suggest  that  they 
may  be  modeled  as  horizontally  propagating  internal  wave  modes  which  are 
near  their  northern  turning  latitude.  In  such  a  model  the  vertical 
structure  of  each  mode  is  independent  of  latitude,  so  the  near  inertial 
vertical  wavenumber  spectrum  is  determined  by  the  high  frequency 
(relative  to  f)  wavenumber  spectrum  at  lower  latitudes.  Fu  (1980), 
following  Munk  and  Phillips  (1968),  shows  that  this  approach  accurately 
models  the  observed  near  inertial  energy  spectra  in  the  deep  North 
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Atlantic.  One  result  of  his  analysis  is  that  the  near  inertial  vertical 
wavenumber  spectrum  should  be  similar  to  that  at  higher  frequencies,  with, 
however,  less  energy  in  the  low  mode  waves,  as  these  have  a  more  south¬ 
erly  turning  latitude.  The  inverse  analysis  of  the  IWEX  experiment 
(Muller  et  al.,  1978),  consistent  with  this  prediction,  finds  little 
energy  at  the  lowest  3  modes  for  near  inertial  motions;  these  modes  con¬ 
tain  significant  energy  only  at  higher  frequencies. 

The  vertical  near  inertial  wave  scales  will  be  estimated  using 
the  GM75  (Garrett  and  Munk,  1975)  wavenumber  spectrum  with  the  addition 
of  a  low  wavenumber  cutoff.  The  wavenumber  spectrum  is  given  by 

A  ( A)  =  - — -  X  >  d  (1) 

[1  +  (X  -  d)]2'5 
=  0  X  <  d 


where  X  =  6/8*  for  vertical  wavenumber  3  and,  6*  =  j*trN/bNQ, 


1.3  km,  N  =  3  cph,  j*  =  6,  and  d  =  j  /j1 
o  p 


The  low  mode  cutoff 


i  is  taken  as  3,  following  Muller  et  al.  (1978).  Using  a  local 

p  - 

Vaisala  frequency  N  =  .45  cph,  the  characteristic  vertical  wavenumber 
6*  l<;  found  to  be  2.2  x  10  1  or  (460  m)  1  .  The  smallest  energetic 

wavenumber  is  j  /j*  $*  or  (920  m)  1 .  Integrating  (1)  45  per  cent  of 
the  energy  is  found  at  wavenumbers  less  than  3*,  90  per  cent  of  the 
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enerqy  is  found  aw  wavenumbers  less  than  43*  =  (115  m)  l.  For  typical 
mixed  layer  height  of  H  =  20  m,  3*H  has  a  value  of  .04,  with  00  per 
cent  of  the  wave  enerqy  having  3»  <  .17.  The  boundary  layer  is  thin  com¬ 
pared  to  tlie  vertical  scale  of  the  energetic  waves. 

Horizontally ,  Fu  (1980)  finds  near-inertial  wave  scales  larger  than 
20  km  for  the  enerqy  containing  waves.  The  interaction  of  the  waves  and 
boundary  layer  will  also  be  modulated  by  the  horizontal  boundary  layer 
scales.  Arm  and  o'Asaro  (1980)  find  gross  boundary  layer  characteristics 
varying  on  scales  10  km  or  larger.  Both  these  scales  are  much  larger  than 
the  boundary  layer  height. 

Equations  of  motion 

The  horizontal  x  momentum  equation  for  an  f  plane,  turbulent, 
Bonus i nesq  fluid  with  no  mean  velocity  can  be  written 
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where  turbulent  stresses  have  been  scaled  by  values  typical  of  a  two- 
dimensional  turbulent  boundary  layer  (Townsend,  1976,  pg .  290)  and  the 
mean  density  has  been  absorbed  into  P.  For  near  inertial  waves  inter¬ 
acting  with  a  turbulent  boundary  layer  of  typical  height  <H  =  20  m)  and 
typical  wave  velocity  (u  —  1  cm/s)  >  length  (L  =  20  km) ,  and  time 
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(T  =  1/f  = 


u 

fL 


.005, 


written : 


4 

10  s)  scales  the  following  non-dimensional  numbers  are  small: 
10  H/L  =  .01.  The  horizontal  momentum  equations  can  now  be 
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The  magnitude  of  P  can  be  estimated  by  Fourier  transforming 
(3)  and  (4)  w  =  /w  e*“ +  ^  and  solving  for  the  kinetic  energy 


(o2  -  f2)2  [  I  u| 2  +  I  vl  2  ]  =  (a2  +  f2)  (k2  +  £2)  I  f> 


2  .  *2  .  M  2  ,  o2i 
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so  that 


uT.  (g2  -  f2) 

(o2  +  f2)*5 


(6) 


Note  that  the  pressure  gradients  above  the  boundary  layer  become  very 
small  near  the  inertial  frequency. 

Armi  and  D'Asaro  (1980)  conclude  that  on  the  time  scale  of  several 
days  or  shorter  the  temperature  fluctuations  observed  at  a  fixed  mooring 
are  primarily  due  to  horizontal  and  vertical  advection.  Turbulent  mass 
transport  becomes  important  only  on  a  longer  time  scale.  Making  this 
assumption  the  appropriate  mass  conservation  equation  for  the  super- 


inertial  motion  is 
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where  f,  is  the  vertical  displacement  and  p  is  the  mean  density.  The 
non-linear  terms  have  been  eliminated  as  from  (2)  above. 

The  vertical  momentum  equation  can  be  written 
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where  N  is  the  Vaisala  frequency,  (7)  has  been  used,  and  the  non-linear 
and  horizontal  stress  divergence  terms  have  been  eliminated  as  in  (2) 
above . 

In  the  absence  of  turbulent  stresses,  equations  (3) ,  (4) ,  (7) ,  (8) , 

and  (9)  are  the  equations  for  inviscid  linear  internal  waves.  For  a 

i$z 

constant  N  they  have  solutions  of  the  form  e  ,  where  p  is  the 
vertical  wavenumber.  In  the  presence  of  stresses  the  pressure  change 
across  the  boundary  layer  is  given  by  the  sum  of  the  two  terms  in  (9) . 

The  vertical  displacement  E,  is  strongly  constrained  by  the  ocean  bot¬ 
tom  at  z  =  0  and  the  internal  wave  displacements  above  the  mixed  layer. 
The  fractional  pressure  change  across  the  boundary  layer  of  thickness 
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H  due  to  the  first  term  in  (9)  should  thus  be  similar  to  that  in  the 
inviscid  case,  BH,  and  thus  be  small.  The  fractional  change  in  pres¬ 
sure  due  to  the  second  bracketed  term  in  (9)  has  a  magnitude 
u*2/ (0  -  f)uL  <  10  using  (6)  to  estimate  P,  and  O  >  1.03f, 
vi*  <  -2  cm/s.  The  fractional  change  in  pressure  across  the  boundary 
ayer  is  thus  small.  This  results  primarily  from  the  large  vertical 
scale  of  the  near  inertial  waves,  relative  to  the  boundary  layer 
thickness. 

Scaled  equations 

The  scaled  equations  for  motions  near  the  inertial  frequency 


are 
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where  the  caret  denotes  a  Fourier  transform  in  time  only.  These 
equations  are  appropriate  for  a  thin  domain  stretching  from  the  ocean 
bottom  z  =  0,  to  just  above  the  boundary  layer.  Note  that  the  stress 
terms  in  (11)  and  (12)  are  the  components  of  turbulent  stress  divergence 
with  frequency  O,  not  the  total  stress  divergences. 


59 


These  equations  can  be  used  to  evaluate  the  interaction  of  the 
boundary  layer  and  the  wave  motions  at  frequency  a.  Consider  a  level 
above  the  boundary  layer  chosen  so  that  the  stress  terms  in  (11)  and 
(12)  are  negligible.  Since  the  pressure  is  constant  across  the  boundary 
layer,  the  pressure  terms  in  (11)  and  (12)  can  be  evaluated  from  the 
velocity  field  at  zq ,  yielding  equations  for  the  boundary  layer 
stresses  in  terms  of  the  velocity  field. 
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where  u  ,  v  are  velocity  components  measured  at  z  . 

o  o  o 

Equations  (15)  and  (16)  state  that  the  presence  of  boundary  layer 
stresses  leads  to  different  velocities  within  the  boundary  layer  than 
above  it,  even  though  the  pressure  field  is  the  same.  They  state  that 
all  observed  velocity  shears  within  the  boundary  layer  are  due  to  tur¬ 
bulent  stresses.  This  is  the  usual  assumption  in  boundary  layer 
analysis.  For  low  frequency  motions,  0  <<  f,  (15)  and  (16)  become 
the  steady  Ekman  layer  equations  with  a  three  way  balance  of  Coriolis 
force,  fu;  pressure  gradient,  f  u^;  and  turbulent  stresses 
(Priestly,  1959,  eg .  34  ).  For  internal  wave  motions,  however,  O  s  f, 
so  the  additional  terms  in  (15)  and  (16)  must  be  included. 
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n*  i  iv  t  •  :  uu  t  .  i 

E;uiti  ns  1 1  !*V  and  (It  I  allow  the  turbulent  stresses  which  act 
;;x>r.  the  motions  at  f  requency  •'  to  be  computed  from  the  velocity  field, 
it  is  precisely  these  dresses  that  may  absorb  energy  from  the  internal 
wavt  motions,  and  t bus  act  to  dissipate  the  waves.  Multiplying  (11), 

111)  and  the  Fourier  transform  of  (9)  by  u  ,  v  and  w  respec¬ 
tively,  (the  daqgei  denotes  complex  comuqat  lun)  ,  summing  the  products, 
aiding  the  con.pl  ex  annulate  of  the  sum  and  using  (14)  ,  an  energy  equa¬ 
tion  is  derived: 

Re  *  7  •  f  i  I  )  -  Re  u  —  I  {T*w"* )  +  \>  ~  ( V  ’w  ’ )  +  w  r~  lw':  )  '  .  ( 1 7 ) 


wh»te  Re  • x:  is  the  real  part  of  x.  Equation  (]7)  states  that  for 
motions  at  frequency  7  the  energy  flux  divergence  is  equal  to  the  rate 
of  working  by  the  boundary  layer  stresses  on  the  velocity  field. 

In  the  Appendix  the  work  done  by  the  third  term  on  the  RHS  of  (17) 
is  shown  to  be  small  compand  to  that  done  by  the  first  two  terms.  The 
turbulent  stress  divergences  can  be  evaluated  usinq  (15)  and  (16)  which 
y  wilds 
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where  C(J)  is  a  complex  velocity  vector,  corresponding  to  clockwise 
motion  for  [Positive  a  and  anticlockwise  motion  for  negative  O. 

C  (j)  is  evaluated  at  the  reference  height  Zq ,  above  the  boundary 
layer . 
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Chapter  8. 


Interpretation  of  near  inertial  observations 


General  features 

The  observed  near  inertial  structure  is  shown  in  figure  11,  Near 
inertial  motions  have  less  energy  within  the  mixed  layer  than  in  the 
interior.  Equations  (15)  and  (16)  require  that  all  differences  in 
velocity  be  due  to  turbulent  stresses.  The  observations  thus  require 
turbulent  stresses  to  fill  the  mixed  layer,  assuming  no  turbulent  stress¬ 
es  in  the  interior.  Note  that  a  random,  statistically  steady  stress 
field  will  have  an  average  inertial  component  of  zero  and  thus  will  not 
affect  the  near  inertial  energy  level.  At  any  time  such  a  stress  field 
is  us  likely  to  increase  the  local  near  inertial  amplitude  as  decrease  it. 
Consistently  less  near  inertial  energy  within  the  mixed  layer  indicates 
a  variation  of  the  boundary  layer  stresses  phase  locked  with  the  near 
inertial  motions,  i.e.,  the  boundary  layer  stresses  are  modulated  by  the 
near  inertial  wave  motions. 

Equation  (17)  gives  the  energy  absorbed  by  the  boundary  layer  from 
the  near  inertial  motions.  Note  that  the  presence  of  stresses  alone 
does  not  lead  to  energy  absorption.  For  this  the  stresses  must  be  cor¬ 
related  with  the  motion  so  that  work  is  done.  Equation  (18)  expresses 
this  same  energy  absorption  as  the  sum  of  terms  involving  the  clockwise 
and  anticlockwise  components  of  the  velocity  field  at  frequency  0. 

Each  term  is  proportional  to  the  quadrature  spectrum  of  the  motion 
within  the  boundary  layer  with  that  at  some  reference  level  outside  the 


boundary  layer.  The  quadrature  spectrum  can  be  expressed  as 


(20) 


where  y  and  0  are  the  coherence  and  phase  of  the  boundary  layer  ve¬ 
locity  with  respect  to  the  reference  velocity.  Since  0  <  y  <  1  the 
sign  of  the  energy  transfer  between  the  wave  velocities  and  the  tur¬ 
bulent  stresses  is  determined  by  the  sign  of  the  phase,  as  shown  in 
table  2.  These  results,  for  the  case  of  energy  absorption  alone,  have 
been  derived  by  Munk  e_t  crK  (1970)  for  a  laminar  Ekman  layer.  Their 
results  are  here  shown  to  be  determined  on  purely  energetic  grounds. 

The  observations  in  figure  11  show  the  boundary  layer  near  iner¬ 
tial  clockwise  velocities  leading  the  interior.  By  table  2  this  implies 
energy  absorption  by  the  boundary  layer  from  internal  wave  motions  with 
0  >  f,  and/or  energy  flux  into  motions  with  0  <  f.  Note  that  although 
the  turbulent  stresses  fill  the  mixed  layer,  they  only  exchange  energy 
with  the  waves  in  the  lower  part  of  the  mixed  layer. 

Energy  flux  calculation 

The  total  work  done  by  the  boundary  layer  stresses  on  the  motions 
at  frequency  J ,  per  unit  area  per  unit  time,  H(C),  is  found  by  inte¬ 
grating  (IB)  across  the  boundary  layer. 
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Table  2 


Energy  flux  and  boundary  layer  phase 
Clockwise  motion 


0  >  f 


B.L.  leads  -  absorption 
B.  l.  lags  "*■  production 


a 

B.L.  leads 
B.L.  lags 


Anticlockwise  motion 


all  0 


B.L.  leads  -*•  absorption 
B.L.  lags  ■*  production 


Lead  and 


lag  are  relative  to  the  direction  of  rotation 


f 

production 

absorption 
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Equation  (21)  is  used  below  to  compute  the  flux  of  near  inertial  wave 
energy  into  the  benthic  boundary  layer. 

For  near  inertial  motions  the  spectral  estimates  C(0,z)  in  (21) 
must  be  evaluated  with  a  high  frequency  resolution  to  resolve  o  -  f. 

C (a )  is  therefore  computed  for  each  instrument  using  a  1024  element 
Fourier  transform  acting  on  a  single  92  day  Hamming  windowed  piece 
(Harris,  197H) .  The  raw  data  is  averaged  and  subsampled  over  127.5 
minutes  before  Fourier  transforming,  resulting  in  a  variance  reduction 
of  approximately  5  per  cent  at  .05  cph.  This  processing  yields  estimates 
of  C(o )  with  a  frequency  resolution  of  approximately  .Olf.  For  a 
given  instrument  pair  Im{  (o  -  *-s  comPute<^  for  each  spectral 

estimate,  these  products  are  then  averaged  over  frequency  to  increase 
ttie  statistical  reliability  of  the  estimate.  This  processing  yields 
spectral  averages  at  a  fixed  height  off  the  bottom  and  ignores  any  vari¬ 
ations  in  mixed  layer  height.  For  a  record  of  only  92  days  length,  this 
is  unavoid. ible.  The  calculation  also  ignores  any  Doppler  shifting  of 
the  wave  frequency  by  the  mean  current. 

The  integrands  in  equation  (21)  are  functions  of  both  frequency 
and  the  pair  of  instruments  used.  The  frequency  dependence  will  be  ex¬ 
amined  first.  In  figure  12b  the  reference  level  z  is  chosen  as  55  m 

o 

and  the  energy  flux  resulting  from  the  clockwise  term  in  (21)  is  plotted 
for  z  =  15  m  as  a  function  of  frequency.  The  solid  line  shows  the 
results;  the  dashed  lines  give  95%  confidence  limits  computed  ^.sing  the 
phase  confidence  limits  given  by  Koopmans  (1974,  pg.  285).  The  computed 


energy  flux  is  indistinguishable  from  zero  at  frequencies  above  .08  cph. 
Some  energy  absorption  is  indicated  at  the  tidal  frequency.  Near 

the  inertial  frequency  significant  energy  absorption  is  indicated  for 
o  >  f;  significant  production  of  energy  is  indicated  for  O  <  f.  Fre¬ 
quencies  below  .02  cph  are  not  well  behaved;  the  magnitude  of  the 
energy  transfer  varies  widely  with  frequency  and  details  of  the  cal¬ 
culation  procedure.  They  are  not  shown  in  figure  12. 

Figure  12a  shows  the  same  computation  with  z  =  65  m,  z  =  55  m. 

o 

Both  these  levels  are  above  the  boundary  layer  and  as  expected  there  is 
no  energy  absorption  indicated.  The  same  calculation  for  the  anticlock¬ 
wise  term  in  (21)  shows  a  weak  pattern  of  energy  absorption  within  the 
boundary  layer  and  near  the  inertial  frequency. 

The  integration  required  by  equation  (21)  is  shown  in  figure  13 
for  two  near  inertial  frequency  bands.  The  results  are  tabulated  in 
table  3.  Figure  13a  shows  the  integration  of  the  clockwise  term  for  a 
band  extending  from  the  inertial  frequency  .04  cph,  to  .069  cph.  Each 
curve  in  13a  plots  the  integrand  from  (21),  i.e.,  the  energy  flux  per 
unit  volume,  as  a  function  of  z  for  a  fixed  reference  height  zq.  The 
area  under  each  curve  gives  the  energy  flux  per  unit  area.  The  solid 
vertical  lines  are  selected  95  per  cent  confidence  limits.  Figure  13a 
indicates  significant  energy  absorption  for  z  <  20  m  for  all  reference 
levels.  The  double  diagonal  line,  a  rough  average  over  all  reference 
levels,  yields  an  energy  absorption  rate  of  .015  ergs/cm  /s.  Notice 
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frequency  cph 


Figure  12 

Knerqy  absorption  density  by  turbulent  boundary  layer  as  computed  from 
equation  (2J)  for  the  clockwise  velocity  component.  The  95%  confidence 
limits  computed  from  phase  confidence  limits  are  shown  as  dashed  lines. 

a)  Reference  instrument  at  65  m,  z  =  55  m; 


b)  Reference  instrument  at  55  m.  z  =  15  m. 
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Figure  13 

Computation  of  near  inertial  energy  absorption  rate  per  unit  area  by 
boundary  layer  using  equation  (21) .  Each  figure  plots  an  integrand  of 
(21)  versus  heiqht  for  5  different  reference  levels,  35  m  to  85  n\.  The 
area  under  each  curve  gives  the  estimated  energy  absorption.  The  clock¬ 
wise  (figures  a,  c)  and  anticlockwise  (figures  b,  d)  contributions  are 
shown  for  two  near  inertial  frequency  bands. 
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however,  that  a  large  contribution  to  this  integral  occurs  for 
z  ■  15  m,  where  there  is  no  data.  Fortunately,  the  contribution  from 

this  region  can  be  bounded  by  noting  that  the  vertical  coordinate  in 
figure  l la  is  a  quadrature  spectrum  and  by  equation  (20)  is  bounded  above 
by  J  (o  -  f ) C (o , z)  ("  (a , zp) |  .  Estimating  this  by  (a  -  f ) | C (o , ) ] 2 , 
z^  15  in,  yields  a  maximum  value  of  the  vertical  coordinate  of 
.018  ergs/cm2 /s  per  10  m.  If  the  diagonal  line  in  13a  is  bent  at 
z^  =  15  so  that  it  intersects  18  at  z  =  0  an  estimate  of 
.024  ergs/cirT/s  for  the  energy  absorption  is  found.  This  value  is 
probably  an  overestimate.  A  lower  bound  can  be  estimated  from  the  area 
under  a  diagonal  line  from  30  m  on  the  horizontal  axis  to 
.003  erqs/cm2/s  per  10  m  on  the  vertical  axis.  This  yields 
.004  ergs/cm2/s. 

Figure  13b  shows  the  same  calculation  for  the  anticlockwise  term 
in  (21).  The  individual  curves  show  a  strong  peak  at  45  m.  As  dis¬ 
cussed  above  the  instrument  at  45  m  is  believed  to  exhibit  a  directional 
bias  of  up  to  10°  relative  to  the  instruments  at  55  m  and  65  m  for  un¬ 
known  reasons.  The  peak  at  45  m  is  likely  due  to  instrumental  effects, 
and  is  not  included  in  the  integration.  This  effect  is  more  pronounced 
for  the  anticlockwise  motions  than  for  the  clockwise  motions,  due  to 
the  smaller  anticlockwise  phases. 


Figures  13c,  d  show  the  same  calculation  for  a  narrower  frequency 
band,  .04-. 049  cph.  The  clockwise  terms,  figure  13c,  have  less  error 
and  larger  coherences  and  phases  than  for  the  wider  band,  figure  13a. 

The  anticlockwise  terms,  figure  13d,  however,  are  barely  distinguish¬ 
able  from  zero  and  show  much  less  order  than  in  figure  13b. 

Internal  wave  energy  absorption 

The  calculation  outlined  above  yields  the  energy  absorption  rates 
shown  in  table  3.  The  results  are  shown  for  two  near  inertial  frequency 
bands,  along  with  estimated  errors.  These  errors  are  dominated  by  the 
lack  of  measurements  below  15  m,  although  the  finite  record  length  also 
makes  a  contribution. 

The  array  used  in  this  experiment  is  short  compared  to  typical 
vertical  internal  vave  scales.  Given  the  instrumental  errors  and  the 
significant  phase  differences  attributed  to  turbulent  stresses,  the 
array  yields  no  information  on  the  vertical  structure  of  the  deep  ocean 
internal  wave  field.  The  spectral  shape  and  energy  level  of  the  ob¬ 
served  superii'.ertial  fluctuations  are  similar  to  the  universal  internal 
wave  forms.  This  alone  suggests  that  these  fluctuations  should  be  at¬ 
tributed  to  internal  waves,  and  the  calculated  energy  flux  be  included 
in  the  internal  wave  field  energy  budget. 

The  above  calculation  yields  energy  fluxes  between  the  boundary 
layer  and  superinertial  and  subinertial,  clockwise  and  anticlockwise 
motions.  Which  of  these  fluxes  should  be  included  in  the  energy  budget 
for  internal  waves?  The  superinertial  clockwise  motions  are  most 
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likely  due  to  internal  waves  and  the  flux  of  energy  associated  with  them 

will  be  included  in  the  internal  wave  energy  budget.  As  noted  in 

chapter  5,  the  near  inertial  anticlockwise  motions  show  far  more  energy 
than  would  be  predicted  by  internal  wave  theory.  This  suggests  that 
these  motions  are  not  due  to  internal  waves.  The  associated  energy  flux 
will  not  be  included  in  the  internal  wave  energy  budget. 

The  above  calculation  suggests  an  energy  flux  into  the  subinertial 
motions,  as  shown  in  figure  12b.  The  change  in  the  sign  of  the  computed 

energy  flux  at  the  inertial  frequency  is  due  to  the  change  in  sign  of 

(0  -  f)  in  equation  (21) ;  the  sign  of  the  relative  phase  of  the  boundary 
layer  and  interior  motions  remains  unchanged.  Motions  with  0  <<  f  are 
believed  to  be  highly  nonlinear  and  quasigeostrophic  (Rhines,  1979); 
motions  with  0  >  f  are  believed  to  be  linear  internal  waves.  The 
dynamics  of  slightly  subinertial  motions  is  unclear.  Since  the  calcu¬ 
lation  used  above  assumes  a  linear  momentum  equation,  it  may  not  cor¬ 
rectly  estimate  the  energy  flux  for  a  <  f. 

The  above  calculation  also  assumes  no  mean  velocity.  Doppler  shif¬ 
ting  of  the  wave  frequency  by  the  mean  velocity  will  result  in  some 
super  inertial  wave  energy  at  subinertial  frequencies.  The  calculation 
will  compute  the  wrong  energy  flux  for  these  waves.  The  computed  energy 
flux  at  slightly  subinertial  frequencies  will  not  be  included  in  the 
internal  wave  energy  budget. 

Only  the  calculated  energy  flux  from  near  inertial  clockwise 
motions  will  be  used  in  the  following  calculations.  This  should  yield 
the  best  estimate  of  the  effect  of  the  benthic  boundary  layer  on  the 
internal  wave  field.  Including  the  anticlockwise  or  subinertial  energy 
absorption  in  the  internal  wave  energy  budget  could  have  a  significant 
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effect  on  the  estimated  energy  fluxes  (c.f.  table  3,  figure  12).  The 


anticlockwise  motions  would  increase  the  total  wave  energy  absorption; 
the  subinertial  motions  would  decrease  it. 

Results  of  calculations 

The  calculated  energy  fluxes  are  shown  in  table  3.  Over  the 
period  of  this  experiment  the  boundary  layer  absorbs  roughly  IT  = 

.015  ergs/cm2/s  from  the  near  inertial  internal  wave  field.  This  is  a 
net  flux.  Although  generation  of  internal  waves  by  the  boundary  layer 
may  be  occurring  (Townsend,  1965) ,  absorption  of  near  inertial  waves 
generated  elsewhere  dominates,  resulting  in  a  net  absorption  of  near 
inertial  energy. 

The  calculated  clockwise  energy  absorption  rate  has  been  used  to 
compute  a  wave  reflection  coefficient 


R(0 )  |  2 


F+(C) 
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n  (o) 
F  (a) 
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where  F  ,  F  are  the  upward  and  downward  internal  wave  energy  fluxes 
respectively.  The  downward  energy  flux  F_  can  be  estimated  from 
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where  C  is  the  vertical  component  of  internal  wave  group  velocity 

g 

(Lighthill,  1978,  chapter  4).  j  R | 2  is  assumed  to  be  close  to  1  so  that 
roughly  half  of  the  internal  wave  energy  E (0 )  is  propagating  downward. 
Equation  (23)  is  evaluated  using  the  same  technique  used  to  evaluate  the 
integrands  of  (21).  A  value  of  3  =  0*  =  (460  m)  roughly  the  value 

of  3  1  weighted  by  the  spectrum  (1),  is  used.  The  resulting  values  of 
F  and  j  R |  2  are  shown  in  table  3.  The  reflection  coefficient  is 
roughly  .95  so  that  5  per  cent  of  the  wave  energy  is  absorbed  in  each 
reflection . 

Another  useful  measure  of  the  energy  absorption  rate  is  the 
boundary  layer  relaxation  time 


_  E(O)  H 

tr  11(0) 


(24) 


which  gives  the  time  required  for  the  observed  energy  absorption  rate  to 
reduce  the  observed  kinetic  energy  within  the  boundary  layer  of  height 
11  to  zero,  if  there  were  no  energy  input.  Unlike  the  reflection  co¬ 
efficient,  T  depends  only  on  boundary  layer  properties  and  does  not 
R 

require  any  knowledge  of  the  wave  properties.  Evaluating  (24)  yields 
relaxation  times  of  several  days. 
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Chapter  9.  clab  model  of  the  boundary  layer 

In  the  previous  section  the  absorption  of  internal  wave  energy  by 
the  benthic  boundary  layer  was  computed.  Below  the  physics  responsible 
for  this  absorption  is  examined.  A  simple  slab  model  with  a  linearized 
quadratic  drag  law  will  be  shown  to  account  for  most  of  the  observed 
internal  wave/boundary  layer  interactions. 

Equations 

Consider  a  mixed  layer  of  height  H,  uniform  in  velocity  and 

density  above  a  flat  bottom,  and  underlying  a  region  with  constant 

Vaisaia  frequency  N  (Figure  14).  A  constant  velocity  U  is  imposed 

in  the  +x  direction.  Above  the  mixed  layer  the  fluid  is  inviscid. 

The  observed  mixed  layers  are  qenerally  capped  by  a  density  change  of 

order  —-1-  ;t  n‘H.  Since  3H  «-<  1  for  vertical  wavenumbers,  typical 

of  the  deep  ocean,  this  density  step  is  dynamically  insignificant 

(P'Asaro,  1978) ,  and  not  included  in  the  model. 

The  observations  presented  here  clearly  indicate  the  presence  of 

turbulent  stresses  throughout  the  mixed  layer.  Their  structure  within 

the  mixed  layer  is  less  clear,  particularly  as  no  measurements  were  taken 

within  1  ‘i  m  of  the  bottom.  An  integral  model  of  the  boundary  layer  will 

thus  be  used.  Within  the  mixed  layer  the  turbulent  stresses,  T,  are 

assumed  to  decay  linearly  from  their  bottom  value  T  to  zero  at  z  =  H. 
‘  o 

;  acts  to  oppose  the  instantaneous  mixed  layer  velocity  and  has  mag- 

o 

ni tilde  ; i  I  •  u*?  =  C  Sr,  where  S  is  the  instantaneous  mixed  layer 

o  D 

speed.  The  linearized  Boussinesq  equations  for  the  perturbation  quan¬ 
tities  in  the  mixed  layer  are: 
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■  lab  model  of  internal  wave  interaction  with  boundary  layer.  The 
layer  is  uni  form  in  density  and  velocity  with  a  drncj  applied  to  it 
the  i -.von  bottom.  The  interior  is  inv.iscid,  with  constant  N,  and 
Townward  rni  mat  inn  wav*'  reflectinq  from  the  mixed  layer . 
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where  r  =  -  results  from  linearizing  the  drag  law,  and  P  denotes 

H 

the  pressure  divided  by  density.  Within  the  mixed  layer  u,  v,  P  and 
3w  3w  w(H) 

- —  are  constant  and  since  w(0)  =  0,  r —  =  - .  The  mixed  layer 

■jZ  d  z  H 

velocities  are  completely  determined  by  the  value  of  P  at  z  =  H. 

This  model  is  similar  to  that  of  Fu  (1980) ,  although  he  does  not  include 
a  mean  current. 

In  the  interior  equations  (25)  -  (27)  apply  with  r  =  0  and 
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These  are  the  well  known  internal  wave  equations  (Garrett  and  Munk, 

.  .  i(kx  +  Sly  -  at) 

1979).  Assuming  that  all  quantities  vary  as  e 

interior  equations  have  a  solution  of  the  form 
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i  -?V*. 


(31) 


p2  =  (k2  +  l2)  (N2  -  O 2 )  /  ( O 2  -  f2) 


This  corresponds  for  (i  >  0  to  a  downward  propagating  internal  wave  of 
unit  amplitude  and  an  upward  propagating  wave  of  amplitude  R.  The 
equations  can  be  solved  by  requiring  P  and  w  to  be  continuous  at 
z  =  H. 

Solving  the  mixed  layer  and  interior  equations  separately  for 
w(H)  in  terms  of  P  yields: 


1  -  R  _ _ a  a2  -  f2  -  2r 2  +  3  iro 

1  +  R  B  H(g2  -  f2)  r(l  +  sin2  0)  -  io 


where  sin  0  =  U,/(k2  +  £2)^. 

Solving  (25)  and  (26)  for  |u|2  +  |v|2  in  terms  of  P  with  and 
without  r,  yields  the  ratio  of  kinetic  energy  within  the  mixed  layer 
to  that  immediately  above  it. 
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(33) 


A  similar  calculation  yields  the  cross  spectra  for  the  clockwise 
component  of  motion 
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(a  -  f)  [a 


-  2rl  +  3ira] 


The  relaxation  time  T  defined  by  (24)  can  be  computed  by  forming  an 

R 

energy  equation  yielding 

_  _  u2  +  f2  -  2rf  sin  8  cos  6  +  r2 (1  +  3  sin7  9) 

R  2r  O2  +  f2  +  a2  cos2  0  +  f2  sin2  0  +  2r2  (1  +  sin2  0) 

Result.:: 

The  slab  model  has  one  parameter,  r  =  U/H.  Using  U/u*  =  32, 

H  =  20  m,  U  =  7  cm/s  yields  r  =  3.4  x  10  6  sec  1 .  Note  that  r  <<  f. 
For  a  given  incident  internal  w.ve  specified  by  a  frequency  0,  a  ver¬ 
tical  wavenumber  (i,  and  a  propagation  direction  0,  the  model  predicts 
the  boundary  layer  velocity  field,  and  a  wave  reflection  coefficient  R. 
Since  both  the  interior  velocity  field  at  z  =  H,  and  the  boundary 
layer  velocity  field  are  determined  by  the  pressure  at  z  =  H,  all 
quantities  involving  just  these  velocities  are  independent  of  the  ver¬ 
tical  wavenumber  of  the  wave. 

Figure  15a  plots  the  relaxation  time  T  (O)  from  equation  (35) . 

It  is  somewhat  larger  than  1  day,  and  roughly  independent  of  frequency 
arid  wave  direction.  This  value  is  somewhat  smaller,  but  certainly  with¬ 
in  the  range  of  values  calculated  from  the  current  meter  data  (Table  3). 


Figure  15 


Slab  model  results  plotted  versus  frequency,  for  wave  orientation 

0  =  -90°,  -45°,  0°,  45°,  90°  and  r  =  3  *  10~6 . 

a)  Relaxation  time  (equation  35) . 

b)  Ratio  of  horizontal  kinetic  energy  density  within  mixed  layer  to 
energy  above  mixed  layer  (equation  33). 

c)  Phase  of  clockwise  velocity  within  mixed  layer  relative  to  interior 
(equation  34). 

d)  Same  for  anticlockwise  velocity. 

e)  Reflection  coefficient  for  wave  energy  assuming  3  =  (460  m)  1 
(equation  32) . 


ACL  Phase  CL  Phase  Energ/  ratio  t  days 


From  equation  (35)  T  is  seen  to  be  roughly  -g—  for  O  a-  f  and 
r/f  <<  1.  This  result  can  be  simply  derived  by  forming  an  energy 
equation  from  equations  (25)  -  (29);  the  energy  absorption  rate  is  seen 

to  be  r  [  2 1  ti  |  +  |v|2).  If  |u|2  =  |  v  [  -  the  energy  absorption  rate  is 

,  ( |  u  |  2  +  |  v  [  2  >  „  ( |  u  |  2  +  |  v  |  2 )  , 

3r  — 1 — 1 - - — 1 — 1 - ,  the  energy  density  is  — 1 — 1 - 1 — 1 so  x  =  l/3r. 

2  2  R 

Figure  15b  plots  the  ratio  of  horizontal  kinetic  energy  in  the 
mixed  layer  with  that  in  the  interior  at  z  =  H  (equation  (33)).  The 
results  are  independent  of  direction,  but  strongly  frequency  dependent. 

For  a  >  1.3f  these  energies  are  nearly  the  same;  for  0  <  1.3f,  the 
boundary  layer  energy  decreases,  reaching  zero  at  0  "  f.  This  behavior 
is  consistent  with  the  observations  which  show  significantly  less  iner¬ 
tial  energy  in  the  bottommost  instruments  only  for  0  near  f. 

A  scale  analysis  of  the  momentum  equations  (25)  and  (26)  explains 
the  special  significance  of  the  inertial  frequency.  Within  the  boundary 
layer  a  balance  of  inertial,  frictional  and  pressure  forces  exists.  The 
magnitude  of  the  pressure  forces  is  (0  -  f)  by  (6),  where  u^  is  the 
typical  velocity  above  the  boundary  layer.  For  o  -  f  >>  r  the  inertial 
forces  have  magnitude  du.  The  frictional  forces,  of  magnitude  ru, 
are  much  smaller  since  r  <<  f.  The  inertial  and  pressure  forces  balance, 
just  as  in  the  interior,  and  the  velocities  inside  the  boundary  layer  are 
the  same  as  above  it.  The  relative  magnitudes  of  the  terms  in  (25)  for 
a  -  f  >>  r  are 
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For  o  -  f  r  this  balance  can  no  longer  hold  and  the  frictional  terms 
become  important.  For  o  -  f  <<  r  they  become  dominant  and  a  new 


balance  exists 


(37) 


The  boundary  layer  velocity  u,  is  much  smaller  than  the  interior 

velocity  u  ,  their  ratio  being  roughly  (0  -  f)/r,  as  in  (33). 

o 

Physically,  the  drag  forces  are  small  and  only  become  important  near  the 
inertial  frequency  where  the  inertial  forces  and  pressure  gradients  in 
tiic  waves  become  small. 

Figures  15c,  d  show  the  relative  phase  of  the  boundary  layer  and 
interior  velocities  for  clockwise  and  anticlockwise  components  respec¬ 
tively.  These  phases  are  small  except  near  the  inertial  frequency, 
for  the  reasons  discussed  above.  Near  the  inertial  frequency  the  boun¬ 
dary  Layer  leads  in  the  clockwise  component  consistent  with  energy  absor¬ 
ption;  the  anticlockwise  phase  is  strongly  dependent  on  0.  The  pre¬ 
dicted  clockwise  phase  ls  roughly  consistent  with  the  observations; 
t  unite  1  ,  for  example,  shows  the  clockwise  phase  leading  for  o  c  1.4f. 
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The  near  inertial  band  anistropy  seen  in  the  data  can  be  interpreted  as 
internal  wave  propagation  parallel  to  the  mean  flow,  i.e.,  9=0.  In 


this  case  the  slab  model  predicts  a  leading  anticlockwise  velocity  com¬ 
ponent  in  the  boundary  layer,  as  observed.  If,  as  seems  likely,  the 
anticlockwise  velocities  are  not  due  primarily  to  internal  wave  motions, 
these  model  predictions  cannot  be  compared  with  the  observations. 

Figure  15e  plots  the  model  reflection  coefficient,  j r| 2 ,  as  a 
function  of  frequency.  Unlike  the  quantities  shown  in  15a  -  d,  | r| 2 , 
depends  on  the  vertical  wavenumber  of  the  incoming  wave.  A  value  of 
0  =  (460  m)  1  is  used,  the  same  value  used  in  Table  3.  The  predicted 
reflection  coefficient  is  significantly  less  than  1  only  near  the 
inertial  frequency.  It  achieves  a  minimum  value  of  approximately  .9, 
at  a  h,  l.lf,  before  rising  to  1  at  o  =  f-  The  value  of  I r| 2  com- 
puted  from  the  data  is  consistent  with  this  behavior,  showing  a  sig¬ 
nificant  absorption  only  near  the  inertial  frequency,  with  a  value  of 
|  r|  2  t  .95  (Table  2,  figure  3).  The  predicted  behavior  of  | r| 2  can  be 
explained  by  noting  that 


H  EML  2 

X  E  C 
R  INT  g 


(39) 


using  (22)  -  (24) .  From  figure  15a  x  can  be  considered  constant; 

R 

E  /E  is  shown  in  figure  15b.  For  a  >>  f  E  /E  =  1.  C  , 

ML  INT  ML  INT  g 

however,  varies  as  o'  -  f2  so  that  1  -  |  R [  2  't  1/ (o2  -  f 2 )  .  For 
j  -  f  <<  r,  E  /E  decreases  as  (a2  -  f2)2  (equation  (33))  so 


1  -  j  Rj  'v,  (a2  -  »  A  minimum  'alue  of  |r|2  therefore  occurs  for 
o  o  f  +  r.  Physically  | rJ 2  is  near  1  for  O  >>  f  because  the 
frictional  terms  are  small.  j  R ) 7  decreases  as  o  approaches  f 
because  the  group  velocity  decreases  and  the  small  frictional  forces 
have  more  time  to  act  on  the  wave.  Very  near  f,  the  frictional 
forces  decrease  the  boundary  layer  velocities  and  (  r|  2  increases 
toward  1 . 

The  least  known  parameter  in  the  above  model  is  8,  the  vertical 
wavenumber.  Noting  that  the  group  velocity  is  proportional  to  6  *, 

(Rj'  is  seen  to  be  highly  dependent  on  8  (equation  39).  Longer  waves 
have  a  faster  group  velocity  and  are  absorbed  less,  while  short  waves, 
conversely,  are  absorbed  more  strongly.  Thus,  although  the  boundary 
layer  may  not  strongly  absorb  the  internal  wave  field  as  a  whole,  its 
effect  may  be  much  stronger  on  the  shorter,  slower  internal  waves. 
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Chapter  10.  Discussion 

The  above  calculation  of  near  inertial  energy  flux  into  the  boun¬ 
dary  layer  is  subject  to  several  uncertainties  not  included  in  the  error 
bounds.  The  array  used  in  this  experiment  is,  as  mentioned  above,  too 
short  to  provide  any  vertical  wavenumber  estimates  for  the  internal  wave 
field.  The  calculation  provides  an  estimate  of  energy  flux  from  the 
near  inertial  motions  into  the  boundary  layer,  but  does  not  require  the 
near  inertial  motions  to  be  related  to  the  internal  wave  field  in  mid¬ 
water.  Leaman  (1976)  finds  a  near  bottom  peak  in  internal  wave  energy 
within  a  few  hundred  meters  of  the  bottom  in  this  area.  This  suggests 
that  the  dynamics  of  the  near  bottom  internal  wave  field  may  be  dif¬ 
ferent  from  that  in  midwater,  so  that  the  flux  measured  here  cannot 
necessarily  be  used  in  the  energy  budget  of  the  midwater  internal  wave 
field.  A  related  source  of  uncertainty  results  from  the  division  of 
the  energy  flux  into  clockwise  and  anticlockwise  sources.  As  discussed 
above  only  the  flux  resulting  from  the  clockwise  near  inertial  motions 
is  used  in  the  calculations.  Similarly,  the  calculation  indicates  a 
significant  energy  flux  out  of  the  boundary  layer  at  slightly  subinertial 
frequencies.  This  flux  has,  perhaps  incorrectly,  not  been  included  in 
the  estimated  near  inertial  wave  energy  flux.  Statistical  estimates  of 
these  uncertainties  are  difficult  to  form,  but  they  probably  will  not 
change  the  sign  or  general  magnitude  of  the  computed  flux. 
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The  estimated  flux  of  energy  from  the  internal  wave  field  given 

above  includes  only  contributions  from  the  near  inertial  motions. 

tnerqy  absorption  from  frequencies  above  .07  cph  has  not  been  included, 

as  the  calculated  absorption  is  below  the  statistical  noise.  The  model 

presented  in  chapter  'a  suggests  that  for  a  »  f  the  absorption  of 

internal  waves  by  the  boundary  layer  can  be  modeled  using  a  constant 

relaxation  time  T  acting  on  the  wave  kinetic  energy.  Using  x  = 

R  R 

2  days,  a  boundary  layer  height  H  =  20  m,  and  an  energy  of  .7  cm2/s2 
for  the  band  .07  -  .45  cph,  an  energy  flux  of  .008  erqs/cm2/s  is 
obtained,  which  adds  significantly  to  the  total  absorption  of  internal 
wave  energy . 

The  relative  importance  of  the  near  inertial  and  low  frequency 

motions,  in  the  boundary  layer  energetics  can  be  estimated  by  comparing 

the  energy  fluxes  from  these  two  sources.  The  low  frequency  turbulent 

ktnetic  energy  production  in  the  boundary  layer  is  roughly  given  by 

integrating  u^w  •  across  the  boundary  layer  with  TT’w"^  u*2  and 
J  z. 

ill  U* 

—  -  —  as  appropriate  for  a  logarithmic  boundary  Layer;  k  .4  is 
von  Karman’s  constant  (Hinze,  1959).  Integrating  from  z  =  4v/u*  to 
a  height  of  20  m  yields  a  total  production  of  .24  ergs/cm'Vs  for 
U  --  7  cm/s.  Thi.s  is  much  more  than  the  energy  absorbed  from  the 
internal  waves.  The  boundary  layer  is  thus  driven  primarily  by  the  low 
frequency  geos trophic  velocities;  the  internal  waves  only  modulate  it. 
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The  flux  of  internal  wave  energy  into  the  boundary  layer  calculated 
here  is  for  one  location  averaged  over  3  months.  The  model  presented  in 
chapter  9  suggests  that  this  flux  should  be  proportional  to  the  near 
bottom  current  speed.  This  predicts  considerably  larger  absorption  of 
the  wave  energy  to  the  north  of  this  site  where  the  currents  are  larger 
(Schmitz,  197G)  and  considerably  less  absorption  for  regions  of  low 
near  bottom  speeds.  The  calculation  is  probably  not  relevant  to  regions 
of  significant  bottom  slope  or  rough  topography  where  the  kinematics  of 
the  internal  wave  field  is  more  complex  (Phillips,  1977,  chapter  5). 

Thus  this  calculation  cannot  be  used  to  estimate  fluxes  in  or  out  of 
the  internal  wave  field  over  rough  topography. 

Previous  observations  of  near  inertial  waves  have  shown  a  net 
downward  energy  flux  (Muller  et  al . ,  1979;  Leaman ,  1976).  if  the  dif¬ 
ference  between  downward  and  upward  going  energy  is  attributed  to  absorp¬ 
tion  m  the  boundary  layer  a  reflection  coefficient  | r] 2  . 5  is 

;reiicted  (Leaman,  1976).  The  analysis  presented  here  is  not  consistent 
with  such  t  small  value  of  |p]‘  and  suggests  that  an  alternate 
expl anat ion ,  such  as  that  advanced  by  Fu  (1980)  is  more  correct. 

The  absorption  of  internal  wave  energy  by  the  boundary  layer  will 
act  as  a  net  sink  of  energy  for  the  internal  wave  field.  Assuming  an 
internal  wave  energy  of  4  <  106  ergs/cm'  (Garrett  and  Munk ,  1979)  a 
flux  of  energy  to  the  boundary  layer  of  .02  ergs  cnT/s  will  absorb  all 
tile  enor-iv  in  file  internal  wave  field  in  roughly  b  years.  For  comparison 
an  internal  wave  driven  vertical  eddy  diffusivitv  of  3  *  10  cm  ,  o 


Bc> 


would  yield  roughly  the  same  wave  energy  absorption  rate  (Garrett  and 
Munk,  1979).  Estimates  of  energy  input  to  the  internal  wave  field  are 
generally  an  order  of  magnitude  larger  than  the  energy  absorption  rate 
calculated  here  (Kase,  1979;  Fu,  1980;  Bell,  1975).  Unless  internal 
wave  field  relaxation  time  is  much  longer  than  indicated  by  these  es¬ 
timates,  the  benthic  boundary  layer  over  a  flat  bottom  plays  a  minor 
i  ill  in  dissipating  internal  wave  energy. 
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